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ABSTRACT 

Tliis thesis provides a comparison of the accuracies of two methods for fitting a 
WeibuU distribution to a set of aircraft engines time-between-failure data. 

One method used is the Maximum Likelihood Method and assumes that these engine 
failure times are independent. The other method is a Modified Method of Moments procedure 
and uses the fact that if time to failure T has a WeibuU distribution with scale parameter X and 
shape parameter P, then has an exponential distribution with scale parameter X^. The latter 
method makes no assumption about independent failure times. 

A comparison is made from times that are randomly generated with a program. The 
program generates times in a manner that resembles the way in which engine failures occur in 
the real world for an engine with three subsystems. These generated operating times between 
failures for the same engine are not statisticaUy independent. This comparison was extended to 
real data. 

Although the two methods gave good fits, the Maximum Likelihood Method produced 
a better fit than the Modified Method of Moments. Explanations for this fact are analyzed and 
presented in the conclusions. 
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THESIS DISCLAIMER 



The reader is cautioned that computer programs developed in this research may not 
have been exercised for all cases of interest. While effort has been made, within the time 
available, to ensure that the programs are free of computational and logic errors, they cannot 
be considered validated. Any application of these programs without additional verification is at 
the risk of the user. 
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EXECUTIVE SUMMARY 



This thesis provides a comparison of the accuracies of two methods for fitting a 
Weibull Distribution to a set of aircraft engines time-between-failure data. 

One method used is the Maximum Likelihood Method and assumes that these engine 
failure times are independent. The other method is a Modified Method of Moments procedure 
and uses the fact that if time to failure T has a Weibull Distribution with scale parameter X and 
shape parameter P, then T*^ has an exponential distribution with scale parameter X^. The latter 
method makes no assumption about independent failure times. 

In a fleet of fielded aircraft engines, engines are repaired when a fault is detected. In 
addition, some critical components are inspected/repaired after a designated number of flight 
hours. When faults are detected the engine is repaired and the component is replaced on the 
engine unless the engine has accumulated a total number of flight hours near to the hard 
inspection time of a critical component(s). If this is the case the critical components(s) are also 
inspected/repaired. In nearly all cases, the engines are not completely overhauled to restore 
them to new condition. After repair, components in the engine have mixed ages with respect to 
flight hours with some components new or repaired. Consequently the times between failures 
on a specific engine are not independent. 

A comparison is made for times that are randomly generated from a program. The 
program generates times in a manner that resembles the way in which engine failures occur in 
the real world for an engine with three subsystems. These operating times between failures for 
the same engine are generated in a way that they are not statistically independent. This is 
accomplished by using a program in S-Plus which simulates failure times for an engine with 
three subsystems. The failure time distributions of each one of these subsystems are Weibull 
wdth fixed parameters determined a-priori and based on real world experience. The simulation 
gathers three failure times for each engine. The study includes failures for two different sets of 
20 and 100 engines respectively. Once this comparison was completed, the study was extended 
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to real data using TF-34 engine operating times between failures from the Naval Aviation 
Logistics Data Analysis (NALDA) data base. 

The estimation of the parameters for the Weibull Distribution which fits the failure 
times was done using spread-sheet software programs. The accuracy of each method was 
measured using the Kolmogorov-Smimov test statistic, embedded in the S-Plus software 
package. 

Although the two methods gave good fits, the Maximum Likelihood Method produced 
a better fit more often than the Modified Method of Moments. Explanations for this fact are 
analyzed and presented in the conclusions. 
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I. INTRODUCTION 



A. BACKGROUND 

1. Engine Reliability Analysis Program (ERAP) 

The Naval Air Warfare Center Aircraft Division Propulsion and Power Engineering 
(NAWCADPPE) was tasked with analyzing aircraft engine failure data to estimate aircraft 
engine reliability and maintainability factors and to assess the effectiveness of existing naval 
aircraft engine maintenance policies. The ERAP program was established to help accomplish 
this task. It generates nine reports designed to identify problem areas and trends in engine 
logistic support, time between engine failures, engine availability, engine readiness degraders, 
and other engine performance factors. The data base used in the ERAP program is the Naval 
Aviation Logistics Data Analysis (NALDA) data base. These reports provide summary 
measures such as mean flight hours between engine failure and other averages and trend lines 
on these averages. They cannot be used, for example, to compute residual reliabilities of 
specified additional lifetimes for engines with known number of repairs and operating times 
since last failure. Nor can they be used to determine optimal maintenance/replacement policies 
for critical engine components. A more detailed analysis of the engine failure data in the 
NALDA data base is required to make these types of assessments. Michael R. Caudill [Ref 1] 
and John A. Malsbury [Ref 2] made more detailed analyses of TF-34 engine performance in 
their Master of Science theses than those provided under the ERAP program. Mark E. Mlikan 
[Ref 3] developed an integral equation for determining optimal no-build times for critical 
components of TF-34 engines relative to their hard inspection times. All three of these theses 
used a Weibull distribution to model the probability distribution of the flight hours between 
engine fidlures for the TF-34 aircraft engines archived in the NALDA data base. 
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2. NALDA Data Base 

NALDA is an operational automated information system. It is a centralized Naval 
Aviation 3M data base supported by the Defense Megacenter, Mechanicsburg, Pa. A user of 
the NALDA data base must make coded queries to the system to obtain specific information. 
This requires considerable experience with NALDA. NALDA contains numerous data 
elements on Navy aircraft engines by engine type and serial number. Anytime an aircraft 
engine is inspected, maintained, removed from an aircraft, repaired, transported, reinstalled or 
discarded, a report is submitted that is translated into the NALDA data base. The NALDA 
data base must be “cleansed” before it can be used to perform some detailed statistical 
analyses such as fitting a probability distribution to the empirical cumulative probability 
distribution of operating times between engine failures. M. Caudill did this cleansing of the 
NALDA data for TF-34 engines provided to him by support staflf at NAWCADPPE. Details 
of this cleansing process and its purpose are described in Caudill’s thesis [Ref 1]. 

3. Operating Times Between Engine Failures 

The NALDA data base reports the operating time between engine removals . Engines 
can be removed for non-failure reasons. Caudill [Ref 1] established a data base of operating 
time between engine failures for TF-34 engines by engine serial number. These times are 
“not” statistically independent because of the manner in which aircraft engines are repaired. 
More often than not, they are not totally overhauled before being reinstalled. Only minor repair 
is made on some occasions. The operating time between the first and second failure of an 
engine is dependent on the operating time to first failure and the extent of the repair. Engines 
with long first failure times and minor first repair will contain many used components with large 
accumulated operating hours when reinstalled. They will likely have shorter operating times to 
the next failure than engines with small operating times to first failure. Consequently, the 
operating times between engine failures, for the same engine, are not statistically independent. 
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B. FITTING THE PROBABILITY DISTRIBUTION OF OPERATING TIME 

BETWEEN FAILURES FOR TF-34 ENGINES 

1. Weibull Distribution 

Aircraft engines wear out. Consequently, any probability distribution used to model 
their time to failure should reflect wear-out characteristics. The Weibull distribution has these 
properties. It was used by Caudill [Ref 1], Malsbury [Ref 2], and Mlikan [Ref 3] to model 
the distribution of operating times between failures for TF-34 engines. They used a modified 
method of moments procedure to compute the parameters of the fitted Weibull distribution. 
This method used special properties of the Weibull distribution and the exponential distribution; 
the sample mean and the sample standard deviation were used to estimate the parameters of the 
fitted Weibull distribution. It was chosen because it was known to be less sensitive to large 
positive or negative correlations (dependency) in the data than the more standard methods of 
maximum likelihood estimation when statistical independence is assumed. 

2. Research Question 

The analysis performed in this thesis provides a comparison of the accuracies of the 
maximum likelihood estimation procedure (assuming independence) and the modified method 
of moments procedure used by Caudill for fitting a Weibull distribution to the probability 
distribution of operating times between engine failures. These comparisons are made for times 
that are randomly generated in a manner that resembles the way in which engine failure times 
occur in the real world. These randomly generated operating times between failures for the 
same engine will not be statistically independent. Once this comparison was completed, it was 
extended to real data using TF-34 engine operating times between failures from the NALDA 
data base. [Ref 1 ] 
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II. GENERATION OF FAILURE TIMES 



A. PROCEDURE 

To generate the failure times, a program in S-Plus (See Appendix A) was developed 
with the following properties. 

1 . Each engine has 3 components (subsystems). 

2. The failure time distributions of each of these components are Weibull distributions 
with fixed parameters determined a-priori and based on experience. 

3. The data generated is a Type II Censored Data Set, and a total of 3 failure times 
are gathered for each engine. 

With these criteria, the program generates the failure times in the following way: 



Step 1: Three Weibull random numbers are generated. They represent the first failure 
time that each one of the three components could have. These values constitute the 
vector faultsl. Three additional Weibull random numbers are generated. They 
represent the second failure time that each one of the three components could have. 
These values constitute the vector _/aM/A2.Three additional Weibull random numbers 
are generated. They represent the third failure time that each one of the three 
components could have. These values constitute the vector faiiltsS. 

Step 2: The minimum value of the vector faultsl is selected, this value is the time of 
the first failure of the engine. The program creates a new vector neufai/2 that is the 
difference between the vector faultsl and the minimum value selected; it is clear that 
this new vector has a component that is 0. This component is replaced in the vector 
ttewfaill by the element of the vector faultsl that has the same index as the first failure. 

Step 3: At this point the engine has two components that have accumulated working 
hours and one that has been just repaired. The vector newfaill has two elements that 
have the remainder of the theoretical first failure for these components and one element 
that is the theoretical second failure for the other component. So, again, the minimum 
element of this vector is selected and it is the second failure time of the engine. 



5 



Step 4: The program creates a new vector called newfailS that is the difference 
between the vector newfail2 and this mirumum value selected; it is clear that this vector 
has a component that is 0. 

If the second engine failure time has taken place in the same component that had the first 
failure, the 0 element is replaced by the element of the vector faultsS whose index matches that 
of the component that has had the two first failures or by the element of the vector faults! 
whose index matches with the component that had the second failure. Again, there is an 
element in the vector faults^ that is a theoretical failure time for one component just repaired, 
and two elements that have the remainder of theoretical f^ures. So, once again the minimum 
element of this vector will be the value of the third engine failure time. The vector enginefails 
gives the three random failure times of the engine. 

B. GENERATION 

An S-Plus program was developed to generate, in the bank of engines, the three failure 
times for the engine. The program can be run for different numbers of engines. In this thesis the 
program was run for 100 engines with a total of 300 failure times and for 20 engines with 60 
failure times. This was done to determine if the number of engines has an apparent effect in the 
results of the study. For each one of these en^ne sets, 30 replications were generated. These 
replications were exported from S-Plus to Excel with the export tool via the S-Plus file-Export 
menu and used to assess averages on the accuracy of the two methods of estimation • 

The first replication for both 20 and 100 engines can be seen in the Appendbc A, all of 
the other replications have been retained on diskettes attached to the master copy of this thesis. 
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HI. ESTIMATION OF PARAMETERS WITH THE MAXIMUM LIKELIHOOD 
AND THE MODIFIED MOMENT METHOD 



A. PROCEDURE 

The two parameter Weibull distribution analyzed in this thesis has a probability density 
function given by: 

/r(0 = . (I) 

or equivalently 

M) = , ( 2 ) 



where ^ ^ ■ Customarily P is called the shape parameter and (or 0 ) is called the scale 



parameter. 

To estimate the parameters, the capabilities of the spread sheet program Excel were 
used. The failure times are generated in S-Plus and exported to Excel, where they are displayed 
automatically on the Sheet 1 Worksheet of Excel. 

The Maximum Likelihood estimators are calculated based on well known formulas 

[Ref 4]. Solve numerically for P ; 



nV Ti^ LnTi 

n " ^ 

- + y LnTi = 0 

P /=i 



1=1 



(3) 



Then; 
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1 




( 4 ) 




(5) 



The parameters for the Modified Method of Moments were calculated based on the 
fact that if T has a Weibull distribution then has a Exponential distribution with a scale 
parameter of . The following demonstration is fi'om [Ref 5]: 

Considering the Weibull distribution given by; 

m = . ( 6 ) 



Let x = t^ so / = , 



dx 




(7) 



Then 



f,(x) = \j\mx )) = 



— X 

P 



1 



1 i-i 

P 



( 8 ) 
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Therefore: 



X ~ exp(/l^ ) . 



(9) 



It is well known that the mean and standard deviation of the exponential distribution 
are equal. Therefore their statistical estimates u and s should be nearly equal. That is 



5 = 



n 



A ; = 1 



1 " 



( 10 ) 



where the 7^ are the times between engine failures. 



The solution for P in equation (10) will be denoted by P . 



1 I 1 " 

Since E{T^) = 6*^ = 1 "T , the sample mean, —V Tf estimates 0 , therefore; 
'•A/ “ 



d = 



1 " 

ntt 



( 11 ) 



This is the same estimate obtained for the MLE of 0. 



B. ESTIMATION 

To carry out the estimation of parameters with the two methods, two Macros in Visual 
Basic were developed. One Macro for 20 engines (See Appendix B), that uses 60 failure times, 
and the other for 100 engines (See Appendix C), that uses 300 failure times. 

Basically, these two Macros are similar and therefore they, and procedures for using 
them, can be similarly described. Once the data is exported to Excel, it appears on Sheet 1. At 
this point it is necessary to open the file with the Macro, called Macro20.xls for the sample of 
20 engines and Macrol00.xls for the sample of 100 engines respectively. Next the window 
with the Macro is hidden from the same Sheet ! the Macro can be run. To run the Macro, 
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choose Macro from the Tools menu. From Macro/Name Reference choose the name of the 
Macro and click on run. The Macro does all the previous calculations and operations needed 
to estimate the parameters and to entitle the entire Worksheet. After these calculations the 
Maero stops. At this point it is neeessary to introduee the following data: 

In Cell B6 introduce the value previously caleulated by the macro for the initial 
value of the shape parameter whieh is in Cell B4. 

This initial parameter was introduced to assure the convergence with the Tool Goal 
Seeker, beeause without it, a solution may not be reaehed. The initial value is given by [Ref 
6 ]; 



Po = 



' 6 

2 


Y^iLnTif-iY^LnTifln 

_ (=1 /=! 




n-\ 




> 



( 12 ) 



Using these choiees for the initial shape parameter, an average of only 3.5 Newton 

iterations were needed to achieve four-place accuracy in the solution for P . 

After introducing the initial value given in Cell B4 into Cell B6, cell C6 must be 
highlighted. Then the Tools Menu is selected, Goal Seek is selected, a 0 is placed in the “To 
Value,” and B6 is entered as the “Changing Cell.” 

Equation (3) is in Cell C6 and the Goal Seeker searches for the value of P which 
satisfies that equation. There were no non-convergence problems in any of the 70 
replications made with that formula. Onee the solution is achieved, the final value for the 
shape parameter appears in Cell B6. The solutions for A, and 0 appear in Cell B9 
(Equation (4)) and B1 1 (Equation (5)) respeetively. These solutions for /3 , X and 0 are 
the maximum likelihood estimates for these parameters of the Weibull distribution that is 
used to model the probability distribution of the family of times between failures on all 
engines. These values are printed on the left side of Sheet2. 

On the right side of Sheet 2 are the ealeulations earried out by the Modified 
Moment Method. 
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There were convergence problems in calculating P using Equation (10) Amazingly 



the same formula may converge or not depending on the terms on each side. Through trial and 
error it was discovered that by writing the formula in the following way it was possible to get 
convergence in all replications: 



n 

1 



In" 

/=1 



, /=1 



(13) 



In Cell Q6 it is necessary to introduce the initial shape parameter given in Cell Q4 
(Equation ( 12 )). Once done, highlight Cell P7 (Equation ( 13 )), select the Tools Menu, select 
Goal Seek and depending on the number of engines, enter in “To Value” 7.811334 (20 
engines, 60 failure times, 60 / ) or 17.349447 (100 engines, 300 failure times, 300 / V 299 ). 
In the “Changing Cell, Q6,” the Goal Seeker looks for the P value that will satisfy Equation 
( 13 ). Cells Q9 and Qll contain the values for the scale parameter X and 0 determined by 
Equation (1 1) and (5). 

A total of 30 replications for each set of engines were done. Afterwards, for the 
reasons explained later, 10 more replications were done using real data. 

Printouts of the failure times and parameter calculations for the 20 and 100 engines first 
replication are provided in Appendix E. 

A summary of the parameters obtained is shown in the following tables: 



11 



Table 1. Twenty Engines = 60 Failure Times 



20 Engines = 60 Failure Times 



Maximum Likelihood Estimator 
Shape Scale 



Modified Moment Method 
Shape Scale 



1 


0.973000000 


251 .705300000 


2 


1.105104390 


279.481900000 


? 


1 .052707000 


242.343600000 


4 


1 .005354000 


282.577600000 


5 


1.211857000 


272.413900000 


0 


1.523761000 


277.443800000 


""7 


1.053815000 


243.377800000 


n 


1.015065000 


330.527000000 


9 


0.961137000 


251 .822700000 


ip 


1.108272000 


263.760200000 




1.258861000 


317.378200000 


12 


1.358459000 


245.953800000 


IP 


1.323305000 


281 .038800000 


14 


1.054978965 


223.276800000 


■A5 


1.228101693 


264.364500000 


tip 


1.161814103 


243.734600000 


17 


1.355617000 


268.067900000 


18 


1.308393000 


309.037200000 


19 


1.391697000 


284.440700000 


20 


1.141031467 


242.289400000 


:2i, 


1 .336786000 


271 132800000 


22 


1.389460000 


244.898600000 




1 .200747000 


288.213200000 


24 


1 .297427000 


293.585500000 


25 


1.400818000 


310.907900000 




1.231510000 


277.128900000 


27 


1.048038281 


277.331100000 




1.093433000 


273.733700000 




1.295133208 


246.789500000 


30 


1.224688831 


303.055100000 



Average 
Std. Dev. 



1 .203679065 272.060400000 

0.14588671 6 25.376732459 



0.985725000 


253.139700000 


1.129925000 


281 .796500000 


1 .036553000 


240.792500000 


1.079900000 


291.030800000 


1.102670000 


263.279100000 


1.442500000 


273.090800000 


1.309674000 


262.946000000 


1.006751000 


329.425200000 


1.078819000 


264.611000000 


1.084296000 


261 .569300000 


1.395249000 


328.224900000 


1.211490000 


237.176300000 


1.461631834 


289.757400000 


1.249979811 


236.965600000 


1.182801000 


260.845800000 


1.239576000 


249.537300000 


1.367488000 


268.798900000 


1 .357252000 


312.609300000 


1 .406588000 


285.359900000 


1.324936416 


255.330700000 


1.223524000 


263.236700000 


1.280865198 


238.618200000 


1.263844000 


293.349300000 


1.260780021 


290.891900000 


1 .325062000 


305 455500000 


1.320218000 


283.439100000 


1 .043436000 


276.830600000 


1 248177000 


286.965800000 


1.240951783 


243.300600000 


1 .402646081 


317.067700000 



1 .235443638 274.848080000 

0. 1 33749260 25.9371 51 577 
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Table 2. One-Hundred Engines = 300 Failure Times 



100 Engines = 

Maximum Likelihood Estimator 
Shape Scale 



300 Failure Times 

Modified Moment Method 
Shape Scale 





1 261028000 


292.556900000 


m- 


1.147177000 


258.382300000 




1.192249000 


267.935700000 




1 191393000 


301 .254500000 


5 


1.100167000 


265 867700000 


6 


1.132621000 


273.303500000 


"r 


1 .220806000 


282.032200000 


8 


1.206135000 


267.017000000 


9 


1 .254734000 


291 .692500000 


10 


1 .284450000 


285.463200000 


11 


1 .277785000 


285.850800000 


"l2 


1 .274025000 


263.121900000 


13 


1 192249000 


267.935700000 


14 


1.194163000 


271 .996400000 


IS 


1.216995000 


272.494100000 


16 


1.169718000 


293.370600000 


17 


1 .283386000 


255.640200000 


M 


1.416971000 


294.799600000 


W 


1 .264662000 


268.682800000 


M 


1.342221000 


279 712400000 


zf 


1 .229806000 


287.109500000 




1.199120000 


267.635600000 


23 


1.189326000 


288.996500000 




1.199541000 


263.550100000 




1 .268509000 


261.778100000 


2 ? 


1.257310000 


256.670600000 


27 


1 .274659000 


270.850700000 


2 ? 


1 208075000 


276 652800000 


29 


1 .220756000 


264.013900000 


M 


1.127942000 


279.203900000 



Average 1.226599300 275.185723333 
Std. Dev. 0.063780678 12.379509005 



1 .288232000 


294 659500000 


1.182793000 


261 .246400000 


1 .250620000 


272.423600000 


1 .243339000 


305 757700000 


1.187671000 


272.5641 00000 


1.205131000 


279 422100000 


1.334481000 


290.504100000 


1 246956000 


270.084700000 


1.271828000 


293.041800000 


1.292318000 


286.039200000 


1 387920000 


293.508300000 


1.341512000 


267.591700000 


1 250620000 


272423600000 


1 .228539000 


274.677200000 


1.231497000 


273.612100000 


1 .200586000 


296.122500000 


1 .259945000 


254 068200000 


1 .422087000 


295.112700000 


1 .379585000 


276.306900000 


1 .448534000 


286 372700000 


1 .269967000 


290.254800000 


1.281756000 


273.804000000 


1.207016000 


290.494900000 


1.346121000 


273.747200000 


1.201983000 


256 990000000 


1 .249969000 


256 167900000 


1.311836000 


273 423600000 


1 .237968000 


278 973700000 


1.219361000 


263 912300000 


1.173691000 


283.295300000 



1.271795400 278.553426667 
0.07051 61 79 1 3.009908323 
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IV. ACCURACY ANALYSIS 



A. PROCEDURE 

The accuracy of each method was measured using the Kolmogorov-Smimov goodness 
of fit Test and the S-Plus software package. The null hypothesis in S-Plus is “the True cdf 
(empirical) equals the Weibull distribution”. This test has the advantage that it uses continuous 
populations. Therefore it is not necessary to group the data. 

This test has the disadvantage that when distribution parameters are estimated from the 
same sample of data, the test produces a conservative result; that is, the probability of a Type I 
error will be smaller than specified [Ref 6]. This drawback has only a relative importance in the 
present study. The K-S Test p-value is a conservative result whether the parameters are 
estimated by the Maximum Likelihood Method or the Modified Moment Method. 

A S-Plus function called Tocompare (See Appendbc E) was created to calculate the K- 
S statistics for the set of parameters calculated in each replication using both methods. This S- 
Plus function computes the K-S Test Statistic {Supremum) and the p-value for each method, 
and it also plots the empirical and hypothesized distributions. In Appendbc E the reader will find 
the output of the function for the first replication of each set of parameters. 

B. ACCURACY 

The calculations for all replications and the values obtained for each set of engines are 
presented in Tables 3 and 4. The letter “ B” marks the method that achieved a better fit, using 
as reference the K-S statistic; that is, the method with the smaller K-S statistic is the more 
accurate procedure for that replication. 

For the set of 20 engines, the MLE method gave a better fit in 17 of the 30 replications. 
For the set of 100 engines, the MLE method gave a better fit in 20 of the 30 replications. 
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Replication Number 



Table 3. Better Fit Twenty Engines 



20 Engines = 60 Failures times 





M.LE. 






M.M.M. 




shape 


scale k-s 


p-value 


shape 


scale k-s 


p-value 



,1 


0.973000000 


251.705300000 


0.054600000 


0.994000000 




0.985725000 


253.139700000 


0.050900000 


0.997700000 


B 


2 


1.105104390 


279.481900000 


0.088300000 


0.737300000 




1.129925000 


281.796500000 


0.081000000 


0.826600000 


B 


•3 


1.052707000 


242.343600000 


0.076300000 


0.876000000 




1.036553000 


240.792500000 


0.073000000 


0.906600000 


B 


4 


1 .005354000 


282.577600000 


0.061700000 


0.976300000 


B 


1.079900000 


291.030800000 


0.086100000 


0.765600000 




5 


1.211857000 


272.413900000 


0.101800000 


0.562800000 




1.102670000 


263.279100000 


0.097300000 


0.621600000 


B 


6 


1.523761000 


277.443800000 


0.086700000 


0.757300000 




1.442500000 


273.090800000 


0.081000000 


0.826400000 


B 


7 


1.053815000 


243.377800000 


0.091500000 


0.696500000 


B 


1 .309674000 


262.946000000 


0.115700000 


0.398400000 




$ 


1.015065000 


330.527000000 


0.052700000 


0.996200000 


B 


1.006751000 


329.425200000 


0.053100000 


0.995900000 




9 


0.961137000 


251.822700000 


0.088000000 


0.742000000 


B 


1.078819000 


264.611000000 


0.098600000 


0.604700000 




to 


1.108272000 


263.760200000 


0.064500000 


0.964400000 




1 .084296000 


261.569300000 


0.061500000 


0.977100000 


B 


tt 


1.258861000 


317.378200000 


0.101500000 


0.566300000 


B 


1.395249000 


328.224900000 


0.126600000 


0.290900000 




^ 12 


1.358459000 


245.953800000 


0.104300000 


0.531000000 


B 


1.211490000 


237.176300000 


0.121500000 


0.338000000 






1.323305000 


281.038800000 


0.075700000 


0.881800000 


B 


1.461631834 


289.757400000 


0.110800000 


0.452600000 




I 14 


1.054978965 


223.276800000 


0.088000000 


0.741900000 




1.249979811 


236.965600000 


0.087200000 


0.751700000 


B 


Z IS 


1.228101693 


264.364500000 


0.070300000 


0.928000000 




1.182801000 


260.845800000 


0.065700000 


0.958200000 


B 


1 16 


1.161814103 


243.734600000 


0.078800000 


0.850400000 


B 


1.239576000 


249.537300000 


0.101500000 


0.567200000 




s 17 


1.355617000 


268.067900000 


0.070200000 


0.929200000 


B 


1.367488000 


268.798900000 


0.072500000 


0.910600000 






1.308393000 


309.037200000 


0.049400000 


0.998600000 


B 


1.357252000 


312.609300000 


0.054400000 


0.994300000 




iS 19 


1.391697000 


284.440700000 


0.083700000 


0.794200000 




1.406588000 


285.359900000 


0.083400000 


0.798500000 


B 


20 


1.141031467 


242.289400000 


0.091800000 


0.692600000 


B 


1.324936416 


255.330700000 


0.116800000 


0.386000000 




21 


1.336786000 


271.132800000 


0.107000000 


0.498400000 




1.223524000 


263.236700000 


0.077500000 


0.863800000 


B 


22 


1 .389460000 


244.898600000 


0.085000000 


0.778500000 


B 


1.280865198 


238.618200000 


0.092700000 


0.681600000 




23 


1.200747000 


288.213200000 


0.083700000 


0.795100000 


B 


1 .263844000 


293.349300000 


0.096600000 


0.629800000 




24 


1.297427000 


293.585500000 


0.058400000 


0.986700000 


B 


1.260780021 


290.891900000 


0.059800000 


0.982700000 




25 


1.400818000 


310.907900000 


0.071100000 


0.921800000 




1.325062000 


305.455500000 


0.071100000 


0.922100000 


B 


26 


1.231510000 


277.128900000 


0.062500000 


0.973200000 




1.320218000 


283.439100000 


0.057400000 


0.989000000 


B 


27 


1.048038281 


277.331100000 


0.078500000 


0.853300000 




1 .043436000 


276.830600000 


0.077100000 


0.867800000 


B 


28 


1.093433000 


273.733700000 


0.078000000 


0.859200000 


B 


1.248177000 


286.965800000 


0.084100000 


0.790300000 




29 


1.295133208 


246.789500000 


0.055600000 


0.992500000 


B 


1.240951783 


243.300600000 


0.067600000 


0.947000000 




30 


1.224688831 


303.055100000 


0.083900000 


0.792200000 


B 


1 .402646081 


317.067700000 


0.095900000 


0.639300000 






mean 


mean 


mean 


mean 




mean 


mean 


mean 


mean 






1.203679065 


272.060400000 


0.078116667 


0.822256667 


17 


1.235443638 


274.848080000 


0.083946667 


0.756066667 


13 




stdev 


stdev 


stdev 


stdev 




stdev 


stdev 


stdev 


stdev 






0.145886716 


25.376732459 


0.015610639 


0.146004782 




0.133749260 


25.937151577 


0.020788358 


0.214131410 
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Replication Number 



Table 4. Better Fit One-Hundred Engines 



100 Engines = 300 Failures times 



M.L.E. 



M.M.M. 





shape 


scale 


k-s 


p-value 




shape 


scale 


k-s 


p-value 




1 


1.261028000 


292.556900000 


0.031900000 


0.920600000 


B 


1.288232000 


294.659500000 


0.038200000 


0.774000000 




2 


1.147177000 


258.382300000 


0.027400000 


0.977600000 


B 


1.182793000 


261.246400000 


0.032700000 


0.905100000 




3 


1.192249000 


267.935700000 


0.047000000 


0.522300000 


B 


1.250620000 


272.423600000 


0.051000000 


0.416900000 




4 


1.191393000 


301.254500000 


0.055800000 


0.306900000 




1.243339000 


305.757700000 


0.046600000 


0.532800000 


B 


5 


1.100167000 


265.867700000 


0.060600000 


0.220000000 




1.187671000 


272.564100000 


0.045600000 


0.561900000 


B 


6 


1.132621000 


273.303500000 


0.043200000 


0.630900000 




1.205131000 


279.422100000 


0.042500000 


0.651400000 


B 


7 


1.220806000 


282.032200000 


0.034300000 


0.873100000 


B 


1.334481000 


290.504100000 


0.048900000 


0.469300000 




B 


1.206135000 


267.017000000 


0.037900000 


0.782600000 




1.246956000 


270.084700000 


0.036300000 


0.823000000 


B 


$ 


1.254734000 


291.692500000 


0.021600000 


0.999000000 


B 


1.271828000 


293.041800000 


0.025500000 


0.989900000 




10 


1 .284450000 


285.463200000 


0.022700000 


0.997900000 




1.292318000 


286.039200000 


0.020800000 


0.999500000 


B 


11 


1.277785000 


285.850800000 


0.044600000 


0.589000000 


B 


1.387920000 


293.508300000 


0.058300000 


0.260200000 




12 


1.274025000 


263.121900000 


0.047400000 


0.511700000 


B 


1.341512000 


267.591700000 


0.048900000 


0.469900000 




1$ 


1.192249000 


267.935700000 


0.029400000 


0.957700000 


B 


1.250620000 


272.423600000 


0.035300000 


0.848500000 




14 


1.194163000 


271.996400000 


0.039600000 


0.735300000 


B 


1.228539000 


274.677200000 


0.048600000 


0.477100000 




16 


1.216995000 


272.494100000 


0.031400000 


0.928900000 




1.231497000 


273.612100000 


0.029400000 


0.957500000 


B 


is 


1.169718000 


293.370600000 


0.034800000 


0.861700000 


B 


1.200586000 


296.122500000 


0.041300000 


0.685300000 




17 


1 .283386000 


255.640200000 


0.031800000 


0.922000000 


B 


1.259945000 


254.068200000 


0.033000000 


0.898900000 






1.416971000 


294.799600000 


0.039800000 


0.729300000 




1.422087000 


295.112700000 


0.038700000 


0.759500000 


B 


Id 


1.264662000 


268.682800000 


0.053800000 


0.349300000 


B 


1.379585000 


276.306900000 


0.065700000 


0.150500000 




20 


1.342221000 


279.712400000 


0.039700000 


0.730800000 


B 


1 .448534000 


286.372700000 


0.067500000 


0.129900000 




21 


1.229806000 


287.109500000 


0.041300000 


0.685900000 




1 .269967000 


290.254800000 


0.040300000 


0.713200000 


B 


n 


1.199120000 


267.635600000 


0.036000000 


0.830800000 


B 


1.281756000 


273.804000000 


0.044400000 


0.596300000 




23 


1.189326000 


288.996500000 


0.032800000 


0.903200000 


B 


1.207016000 


290.494900000 


0.033600000 


0.886800000 




24 


1.199541000 


263.550100000 


0.032000000 


0.919200000 


B 


1.346121000 


273.747200000 


0.048900000 


0.469200000 




25 


1.268509000 


261.778100000 


0.028200000 


0.971000000 


B 


1.201983000 


256.990000000 


0.041500000 


0.680900000 




26 


1.257310000 


256.670600000 


0.026400000 


0.984800000 


B 


1.249969000 


256.167900000 


0.027700000 


0.975800000 




27 


1 .274659000 


270.850700000 


0.038700000 


0.759500000 


B 


1.311836000 


273.423600000 


0.048400000 


0.482200000 




23 


1.208075000 


276.652800000 


0.035200000 


0.851300000 




1.237968000 


278.973700000 


0.028200000 


0.971300000 


B 


29 


1.220756000 


264.013900000 


0.033800000 


0.882800000 


B 


1.219361000 


263.912300000 


0.034200000 


0.874100000 




30 


1.127942000 


279.203900000 


0.035700000 


0.838800000 




1.173691000 


283.295300000 


0.031400000 


0.928300000 


B 




mean 


mean 


mean 


mean 




mean 


mean 


mean 


mean 






1 .226599300 


275.185723333 


0.037160000 


0.772463333 


20 


1.271795400 


278.553426667 


0.041113333 


0.677973333 


10 




standev 


standev 


standev 


standev 




standev 


standev 


standev 


standev 






0 063780678 


12.379509005 


0.009086679 


0.209187169 




0.070516179 


13.009908323 


0.010947594 


0.247541499 
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If the MLE and MMM methods are equally good then one would expect that the 
probability the MLE gives a better fit ( using the K-S statistic) should be 0.5. Let represent 
the probability the MLE gives a better fit; the above sets of data can then be used to test Ho : n 
= 0.5 versus Ha ; k ^0.5. For the set of 20 engines, 17 successes were observed out of 30 
replications ( MLE gave a smaller K-S value than MMM); the p- value for testing Ho : 7C = 0.5 
is 0.5847 so Ho would be accepted for any a < 0.5847. For the set of 100 engines, 20 of the 
30 replications gave successes; in this case the p-value is 0.0987 and Ho is rejected for any a 
>0.0987. Certainly this comparison of MLE and MMM, using the binomial test, shows no real 
difference for 20 engines and marginally indicates superiority of MLE for 100 engines. 

If the K-S statistic mean and standard deviation are analyzed it is clear that the MLE. 
gives the smaller mean and standard deviation of the two cases. 

For this data, it seems that despite the fact that there is some dependence between the 
failure times, the MMM does not give a better fit than the MLE. This is true at least for the 
failure times generated with the simulation program used in this study. This is likely due to the 
fact that the failure times for different engines (among engines) are independent. 

At this point real failure data was examined. This real data is fi'om [Ref 1], and fi'om 
the cleansed data collected in Appendix B (Five-Year Time Between Failures Data Base). This 
data was cleansed, removing zeros as failure times. Ten samples of data were taken fi'om the 
data set of operating times between failures on TF 34 aircraft engines. Five samples were of 60 
failures times and five samples were of 300 failure times. 

The same procedure to estimate the parameters and analyze the accuracy of the results 
were applied. The results are shown in Tables 5 and 6. 

With the real data the result is the same, a better fit is obtained more often with the 
MLE method than with the MMM method. In the 60 failure times case the MLE method was 
the better procedure in four of the five samples, and in the 300 failure times case it was better 
three of the five samples using the K-S statistic as measure of accuracy. 
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It is important to point out that in these real data, the engine structure is unknown. 
There is no information about how many elements the engine has, or of course the way that 
these components failed. 

Table 5. Estimate Parameters and Analyze Accuracy Twenty Engines 



20 Engines = 60 Failures times 



M.LE. 



M.M.M. 





shape 


scale 


k-s 


p-value 


shape 


scale 


k-s 


p-value 


1 


1.554927 


651.243400 


0.103000 


0.547400 B 


1.732451 


669.495500 


0.103200 


0.544600 




1.435165 


710.748800 


0.079700 


0.840000 B 


1.735014 


748.435600 


0.130200 


0.261300 


3 


1.531311 


603.432400 


0.072600 


0.909800 


1.635563 


613.677100 


0.051800 


0.997100 B 


4 


1.436419 


609.475100 


0.103800 


0.537100 B 


1.485100 


615.440400 


0.114900 


0.406200 


5 


1.535733 


579.331200 


0.135400 


0.221300 B 


1 .490634 


574.988000 


0.145900 


0.155400 




mean 


mean 


mean 


mean 


mean 


mean 


mean 


mean 




1.414134 


642.431620 


0.048360 


0.511380 4 


1.587754 


663.677300 


0.057660 


0.289300 1 




stdev 


stdev 


stdev 


stdev 


stdev 


stdev 


stdev 


Stdev 




0.045711 


14.890226 


0.008406 


0.185231 


0.023623 


16.834529 


0.006280 


0.138620 



.i: a> 

o E 

g-z. 

a 



C. ANALYSIS OF THE RESULTS 

The author expected “ a-priori” a better fit with the MMM than with the MLE, based 
on the fact that the generated operating times between failures for the same engine are not 
statistically independent. But the resulting fit for the data was better with the MLE more often 
than with the MMM. This result is further analyzed here, by examining how the dependence 
between failure times affect each method. 

To assess dependence between failure times, the correlation and the grade of 
dependence between the failure times was calculated for each replication. The correlation 
coefficients were calculated and added at the end of each replication (See Appendbc A). In 
addition the Kendall-Tau independence test was performed. The null hypothesis in S-Plus for 
this test is “ Tau is equal to 0” that is the two samples are independent. This analysis was 
performed for all possible pairwise combinations of the failure times of an engine; namely first- 
second, first-third, second-third correlations, tau coefficients and p- values were calculated. For 
the Kendall-Tau test the minimum, maximum and mean p-values of those pairwise 
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combinations of each replication were added to the tables. The results are shown in Tables 7 
through 10. 

Table 6. Estimate Parameters and Analyze Accuracy One-Hundred Engines 



100 Engines = 300 Failures times 



M.LE. 



M.M.M. 





shape 


scale 


k-s 


p-value 




shape 


scale 


k-s 


p-value 




1 


1.462681 


618.916600 


0.047200 


0.515100 




1.579829 


632.033000 


0.045700 


0.559000 


B 


2 


1.439461 


663.159800 


0.045200 


0.571700 


B 


1.586168 


681.405000 


0.063000 


0.184300 




3 


1.366409 


641.062700 


0.040500 


0.709200 


B 


1.564422 


666.063600 


0.057300 


0.279000 




4 


1.351920 


636.911100 


0.064600 


0.163100 




1.575530 


665.621100 


0.061900 


0.200100 


B 


5 


1.450201 


652.107900 


0.044300 


0.597800 


B 


1.632821 


673.263800 


0.060400 


0.224100 






mean 


mean 


mean 


mean 




mean 


mean 


mean 


mean 






1.414134 


642.431620 


0.048360 


0.511380 


3 


1.587754 


663.677300 


0.057660 


0.289300 


2 




stdev 


stdev 


stdev 


stdev 




stdev 


stdev 


stdev 


stdev 






0.045711 


14.890226 


0.008406 


0.185231 




0.023623 


16.834529 


0.006280 


0.138620 





.2 5 

o C, 
3 

g-z: 



From those tables may be deduced that with the set of 20 engines the correlation is 
clearly larger than with 100 engines. When the independence is tested with the Kendall-Tau 
test, the values of the tau’s coeflScients are similar to the correlation coefficients, and have 
larger values for the set of 20 engines than the set of 1 00 engines. In general for this set of data, 
the correlation coefficients are small and the p-values in the Kendall-Tau are not significant 
enough to reject the null hypothesis of independence at the usual values. That means the 
dependence obtained in the simulation is small, specially with the set of 100 engines; this may 
be unfavorable to the performance of the MMM. 

To summarize the results, Figures 1 through 4 present the min and max Kendall-Tau 
test p-value for each replication for the method that gives the best fit. 

From these figures and specially from the max p-values, it may be deduced that, for this 
data, in general the larger the max p-value (more independence) the better the fit with the MLE 
method. Conversely the smaller the max p-value (less independence) the better the fit with the 
MMM method. 
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Table 7. Correlation - 20 Engines 





Correlation Coefficients for 20 ennines 






r«p1 


Failurel 


Failure2 


Failures 




rep16 


Failurel 


Failure? 


Failures 






Failurel 


1 








Failurel 


1 










Failure2 


-0 099201099 


1 






Failure2 


0 381725000 


1 








Failures 


-0 325682845 


0 177625785 


1 




Failures 


0 095575000 


-0 033370000 


1 






rep2 


Failurel 


Failure2 


Failures 




rep17 


Failurel 


Failure? 


Failures 






Failurel 


1 








Failurel 


1 










Failure2 


-0 251973035 


1 






Failure2 


-0 214350000 


1 








Failures 


0 460278255 


-0.055862733 


1 




Failures 


-0 567440000 


0 066732000 


1 






ropS 


Failurel 


Failure2 


Failures 




rep18 


Failurel 


Failure? 


Failures 






Failurel 


1 








Failurel 


1 










Feulure2 


-0 292364468 


1 






Failure? 


-0 019240000 


1 








Failures 


-0 171745200 


-0 161919234 






Failures 


0 089556000 


0 366561000 


1 






rep4 


Failurel 


Failure2 


Failures 




rep19 


Failurel 


Failure? 


Failures 






Failurel 


1 








Failurel 


1 










Failure2 


-0 266964793 


1 






Failure? 


0 062676000 


1 








Failures 


0 375627368 


-0 020033442 


1 




Failures 


-0 211740000 


-0 000530000 


1 






repS 


Failurel 


Failure2 


FeulureS 




rep20 


Failurel 


Failure? 


Failures 






Failurel 


1 








Failurel 


1 










F2ulure2 


-0 125786346 


1 






Failure? 


-0 227550000 


1 








Failures 


-0 240386929 


-0 146477812 


1 




Failures 


0 225389000 


-0 371290000 


1 






rep6 


Failurel 


Failure2 


Failures 




rep21 


Failurel 


Failure? 


FulureS 






Failurel 


1 








Failurel 


1 










Failure2 


-0 049045517 


1 






Failure? 


-0 267500000 


1 








Failures 


0 141463053 


0 466711164 


1 




Failures 


0 140151000 


-0 088550000 


1 






rep7 


Failurel 


Failure2 


Failures 




rep?? 


Failurel 


Feulure? 


Failures 






Failurel 


1 








Failurel 


1 










Fzulure2 


-0 334293594 


1 






Failure? 


-0 237640000 


1 








Failures 


-0 102335000 


0 108178021 


1 




FwlureS 


0 051596000 


-0 117030000 


1 






rep8 


Failurel 


Failure2 


Failures 




rep23 


Failurel 


Failure? 


Failures 






Failurel 


1 








Feulurel 


1 










Failure2 


-0 362844249 


1 






Failure? 


-0 300890000 


1 








Failures 


-0 153316474 


-0 014631447 


1 




Failures 


0 077012000 


-0 043760000 


1 






rep9 


Frulurel 


Failure2 


Failures 




rep24 


Failurel 


Failure? 


Failures 






Failurel 


1 








Failurel 


1 










Failure2 


0 138519930 


1 






Failure? 


-0 159170000 


1 








Failures 


-0 383549222 


0 096601004 


1 




Failures 


0 008203000 


0 0873060X 


1 






rep10 


Faulurel 


Failure2 


Failures 




rap?5 


Failurel 


Failure? 


Failures 






Failurel 


1 








Failurel 


1 










Failure2 


-0 205353856 


1 






Failure? 


0 368302000 


1 








Failures 


-0 202946234 


-0 062062633 


1 




Failures 


-0 205760000 


-0 176000000 


1 






repll 


Failurel 


Failure2 


Failures 




rep26 


Failurel 


Failure? 


Failures 






Failurel 


1 








Failurel 


1 










Failure2 


-0 23980229 


1 






Failure? 


-0 076760000 


1 








FeulureS 


0 164227499 


-0 089182145 


1 




Failures 


0 025953000 


-0 17618 


1 






rep12 


Failurel 


Failure2 


Failures 




rep27 


Failurel 


Failure? 


Failures 






Failurel 


1 








Failurel 


1 










Failure2 


-0 291025767 


1 






Failure? 


-0 059610000 


1 








Failures 


-0 283057458 


-0 216673298 


1 




Failures 


0 009859000 


0 421248000 


1 






rep1S 


Failurel 


Failure2 


Failures 




rep28 


Failurel 


Failure? 


Failures 






Failurel 


1 








Failurel 


1 










Failure2 


0 010730000 


1 






Failure? 


-0 166470000 


1 








Failures 


-0 148510000 


0 018339000 


1 




Failures 


-0 408460000 


-0 182890000 


1 






rep14 


Failurel 


Faiture2 


Failures 




rep?9 


Failurel 


Failure? 


Failures 






Failurel 


1 








Failurel 


1 










Failure2 


-0 431010000 


1 






Failure? 


-0 183030000 


1 








Failures 


-0 284860000 


-0 144210000 


1 




Failures 


-0 007070000 


-0 235020000 


1 






rep15 


Failurel 


Faiture2 


Failures 




repSO 


Failurel 


Failure? 


Failures 






Failurel 


1 








Failurel 


1 










Failure2 


-0 075920000 


1 






Failure? 


0 011941000 


1 








Failures 


0 111095000 


-0 381490000 


1 




Failures 


0 144591000 


0 137677000 


1 
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Table 8. Correlation - 100 Engines 



Correlation Coefficients for 100 enoines 



rep1 


Failurel 


Failire2 


FailireS 




rep16 


Failirel 


FailireS 


Failures 


Failirel 


1 








Failirel 


1 






Failire2 


0 012345873 


1 






Failure? 


-0 34S673022 


1 




Fatlir*3 


-0 151985402 


0 027275574 


1 




FailireS 


-0 065777779 


-0 047403929 


1 


rep2 


Failu^el 


Failtre2 


Failures 




rep 17 


Failirel 


FailireS 


Failire3 


Faitirol 


1 








Failirel 


1 






Failira2 


-0 053286726 


1 






Failire2 


-0 095841896 


1 




FailtraS 


-0 155893588 


-0 044274299 


1 




FailireS 


0 153445659 


-0 210881943 


1 


rap3 


Failurel 


Failire2 


FailireS 




rsp18 


Failirel 


FailireS 


Failire3 


Failia-el 


1 








Failirel 


1 






Failure2 


-0 250749723 


1 






Failire2 


0 052473163 


1 




FailiraS 


0 062100274 


-0 094868699 


1 




FailireS 


0 041985512 


0 022318537 


1 


rep4 


Failirel 


Failire2 


FailtreS 




rep19 


Failirel 


FailireS 


FailireS 


Failirel 


1 








Failirel 


1 






Faiiia’«2 


-0 218109532 


1 






Failures 


-0 062164561 


1 




Failire3 


-0 144770967 


-0 109527904 


1 




FailireS 


-0 133619298 


-0 160869201 


1 


repS 


Failirel 


Failire2 


FailireS 




rep20 


Failirel 


Failures 


FailireS 


FailLral 


1 








Failirel 


1 






Faili«^*2 


-0 150008998 


1 






FailireS 


-0 083449627 


1 




FailLT«3 


0 132158154 


-0 078813135 


1 




FailireS 


-0 003128899 


0 055670639 


1 


rep6 


Failirel 


Failire2 


FailireS 




repSI 


Failirel 


FailireS 


FailireS 


Faiiifal 


1 








Failurel 


1 






Fatlir*2 


-0.079147073 


1 






FailireS 


-0 110695458 


1 




Failures 


-0 240225116 


-0 170526078 


1 




FailireS 


-0 170345154 


-0 103486186 


1 


rap7 


Failirel 


Failire2 


FailireS 




rep22 


Failirel 


FailireS 


FailireS 


Failurel 


1 








Failirel 


1 






Failire2 


-0 065957419 


1 






FailireS 


-0 180701469 


1 




FailtreS 


-0 065309916 


-0 149967098 


1 




FailireS 


0 156143924 


-0 219732137 


1 


rap8 


Failurel 


Failire2 


Failures 




rep S3 


Failirel 


FailireS 


FailireS 


Failurel 


1 








Failirel 


1 






Failiae2 


-0 087122634 


1 






FailireS 


-0 073156230 


1 




Failire3 


0 009267376 


-0 198193315 


1 




FailireS 


-0 171758743 


-0 083586000 


1 


rep9 


Failurel 


Failire2 


FailireS 




rep24 


Fiulirel 


Failures 


FailireS 


FailLrel 


1 








Failirel 


1 






Pailu'e2 


-0 024708152 


1 






FailireS 


-0 063590173 


1 




FailixeS 


-0 088238686 


-0 153996216 


1 




FwlireS 


-0 106654719 


-0 162453112 


1 


rep10 


Failirel 


Failire2 


FailireS 




repSS 


Failirel 


FailireS 


Faitire3 


Failirel 


1 








Failirel 


1 






Failu^e2 


-0 085599667 


1 






Failures 


-0 084052612 


1 




FailLT»3 


-0 208242972 


-0 04475226 


1 




Failures 


0 049055891 


-0 063152872 


1 


rep 11 


Failirel 


Failire2 


Failures 




rep26 


Fwlirel 


FailireS 


FailireS 


Failurel 


1 








Failirel 


1 






FailLre2 


-0 030401682 


1 






FailireS 


-0 066404142 


1 




Failtre3 


0 019141802 


0 012598681 


1 




FailireS 


0 099904638 


♦0 219739207 


1 


rep12 


Failirel 


Failire2 


Failures 




rep 27 


Failirel 


Failtre2 


Failire3 


Failiael 


1 








Failirel 


1 






Failu^a2 


-0 243204236 


1 






FailireS 


-0 056026419 


1 




FailiaeJ 


0 016811427 


-0 06918138 


1 




FailireS 


-0.077399648 


-0 002855117 


1 


rep13 


Failirel 


Failire2 


Failures 




repSS 


Failurel 


FailireS 


FwlireS 


Failiael 


1 








Failirel 


1 






Failtre2 


-0 250749723 


1 






FailireS 


-0 081782938 


1 




FailiaeS 


0 062100274 


•0 094868699 


1 




FailireS 


-0 061450025 


0 020642837 


1 


rep14 


Failirel 


Failire? 


Failures 




rep29 


Failurel 


FailireS 


FailireS 


Feiilurel 


1 








Failirel 


1 






FajlLTe2 


-0 064039892 


1 






FailireS 


-0 243170112 


1 




FailtreS 


-0 077445912 


-0 077941408 


1 




FailireS 


0 093441256 


-0 10701961 


1 


repIS 


Failirel 


Failire2 


FailireS 




rep30 


Failirel 


FailireS 


FailireS 


Failtrel 


1 








Failirel 


1 






FailLre2 


-0 098178961 


1 






FailireS 


-0 23939123 


1 




FaiUreS 


0 077241660 


-0 170746677 


1 




FailireS 


-0 020603539 


-0 257213197 


1 
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Table 9. Kendall -Tau Test - 20 Engines 



KENDALL-TAU TEST OF INDEPENDENCE FOR THE SET OF 20 ENGINES 





tau 12 


pvalue 


tau13 


pvalue 


tau23 


pvalue 


mean 


min 


max 


best method 


1 


-0.084 


0.604 


-0.242 


0.136 


0.126 


0.436 


0.392 


0.136 


0.604 


m.m.m 


2 


-0.211 


0.194 


0.221 


0.173 


-0.062 


0.697 


0.355 


0.173 


0.697 


m.m.m 


3 


-0.116 


0.475 


-0.158 


0.330 


-0.095 


0.559 


0.455 


0.330 


0.559 


m.m.m 


4 


-0.179 


0.270 


0.105 


0.516 


0.126 


0.436 


0.408 


0.270 


0.516 


m.I.e. 


5 


-0.126 


0.436 


-0.295 


0.069 


-0.074 


0.650 


0.385 


0.069 


0.650 


m.m.m 


6 


0.105 


0.516 


0.105 


0.516 


0.221 


0.173 


0.402 


0.173 


0.516 


m.m.m 


7 


-0.232 


0.153 


-0.011 


0.948 


0.084 


0.604 


0.568 


0.153 


0.948 


m.I.e. 


a 


-0.221 


0.173 


-0.063 


0.697 


-0.168 


0.299 


0.390 


0.173 


0.697 


m.I.e. 


9 


0.063 


0.697 


-0.295 


0.069 


0.179 


0.270 


0.345 


0.069 


0.697 


m.I.e. 


10 


-0.011 


0.948 


-0.274 


0.092 


-0.126 


0.436 


0.492 


0.092 


0.948 


m.m.m 




-0.147 


0.364 


0.010 


0.948 


-0.063 


0.697 


0.670 


0.364 


0.948 


m.I.e. 


12 


-0.116 


0.475 


-0.158 


0.330 


-0.095 


0.559 


0.455 


0.330 


0.559 


m.I.e. 


.1:;3 


-0.042 


0.795 


-0.179 


0.270 


0.147 


0.364 


0.476 


0.270 


0.795 


m.I.e. 


14 


-0.263 


0.105 


-0.179 


0.270 


-0.179 


0.270 


0.215 


0.105 


0.270 


m.m.m 


'■16 


-0.168 


0.299 


0.232 


0.153 


-0.326 


0.044 


0.166 


0.044 


0.299 


m.m.m 


16 


0.189 


0.243 


0.221 


0.173 


-0.021 


0.897 


0.438 


0.173 


0.897 


m.I.e. 


17 


-0.126 


0.436 


-0.400 


0.014 


0.116 


0.475 


0.308 


0.014 


0.475 


m.I.e. 


18 


0.105 


0.516 


0.053 


0.746 


0.105 


0.746 


0.669 


0.516 


0.746 


m.I.e. 


19 


-0.126 


0.436 


-0.063 


0.697 


0.011 


0.948 


0.694 


0.436 


0.948 


m.m.m 


20 


-0.116 


0.475 


0.221 


0.173 


-0.200 


0.218 


0.289 


0.173 


0.475 


m.I.e. 


21 


-0.168 


0.218 


0.274 


0.092 


-0.095 


0.559 


0.289 


0.092 


0.559 


m.m.m 


22 


-0.211 


0.194 


-0.095 


0.559 


0.042 


0.795 


0.516 


0.194 


0.795 


m.I.e. 


23 


-0.189 


0.243 


0.242 


0.136 


-0.042 


0.795 


0.391 


0.136 


0.795 


m.I.e. 


24 


-0.074 


0.650 


-0.116 


0.475 


0.032 


0.846 


0.657 


0.475 


0.846 


m.I.e. 


25 


0.147 


0.364 


-0.253 


0.119 


-0.053 


0.746 


0.410 


0.119 


0.746 


m.I.e. 


26 


-0.042 


0.795 


0.021 


0.897 


-0.053 


0.746 


0.813 


0.746 


0.897 


m.m.m 


27 


-0.074 


0.650 


0.147 


0.364 


0.105 


0.516 


0.510 


0.364 


0.650 


m.m.m 


28 


-0.042 


0.795 


-0.263 


0.105 


-0.084 


0.604 


0.501 


0.105 


0.795 


m.I.e. 


29 


-0.074 


0.650 


-0.105 


0.516 


-0.105 


0.516 


0.561 


0.516 


0.650 


m.I.e. 


30 


-0.095 


0.559 


0.032 


0.846 


0.032 


0.846 


0.750 


0.559 


0.846 


m.I.e. 




mean 

-0.08804 

stdev 

0.11499 


mean 

0.45767 

stdev 

0.22407 


mean 

-0.04211 

stdev 

0.19334 


mean 

0.38100 

stdev 

0.29522 


mean 

-0.01714 

stdev 

0.12607 


mean 

0.55823 

stdev 

0.23159 


mean 

0.46563 

stdev 

0.15545 


mean mean 
0.24566 0.69414 
stdev stdev 
0.18088 0.18533 
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Table 10. Kendall -Tau Test - 100 Engines 



u 

CJ 

s 

3 

z 

s 

.2 

u 

* 0 . 

X 



KENDALL-TAU TEST OF INDEPENDENCE FOR THE SET OF 100 ENGINES 





tau12 


pvalue 


tau13 


pvalue 


tau23 


pvalue 


mean 


min 


max 


best method 


1 


0.029 


0.668 


-0.136 


0.045 


-0.072 


0.292 


0.335 


0.045 


0.668 


m.I.e. 


2 


-0.093 


0.173 


-0.181 


0.008 


0.003 


0.967 


0.382 


0.008 


0.967 


m.I.e. 


3 


-0.181 


0.008 


0.003 


0.967 


-0.025 


0.708 


0.561 


0.008 


0.967 


m.I.e. 


4 


-0.144 


0.034 


-0.792 


0.428 


-0.064 


0.347 


0.270 


0.034 


0.428 


m.m.m 


5 


-0.089 


0.190 


0.088 


0.192 


-0.036 


0.592 


0.325 


0.190 


0.592 


m.m.m 


6 


-0.138 


0.042 


-0.166 


0.015 


-0.098 


0.150 


0.069 


0.015 


0.150 


m.m.m 


7 


-0.028 


0.681 


-0.059 


0.381 


-0.064 


0.344 


0.469 


0.344 


0.681 


m.I.e. 


6 


-0.036 


0.592 


-0.053 


0.435 


-0.145 


0.033 


0.353 


0.033 


0.592 


m.m.m 


9 


-0.025 


0.708 


-0.009 


0.896 


-0.109 


0.108 


0.570 


0.108 


0.896 


m.I.e. 


10 


-0.026 


0.703 


-0.110 


0.104 


-0.019 


0.775 


0.527 


0.104 


0.775 


m.m.m 


11 


-0.002 


0.976 


0.042 


0.536 


0.026 


0.703 


0.738 


0.536 


0.976 


m.I.e. 


12 


-0.132 


0.052 


0.064 


0.347 


-0.017 


0.798 


0.399 


0.052 


0.798 


m.I.e. 


13 


-0.181 


0.008 


0.003 


0.967 


-0.025 


0.708 


0.561 


0.008 


0.967 


m.I.e. 


14 


-0.081 


0.231 


-0.069 


0.308 


-0.018 


0.789 


0.443 


0.231 


0.789 


m.I.e. 


15 


-0.051 


0.449 


0.095 


0.162 


-0.093 


0.173 


0.261 


0.162 


0.449 


m.m.m 


16 


-0.250 


0.000 


-0.060 


0.378 


-0.039 


0.568 


0.315 


0.000 


0.568 


m.I.e. 


17 


-0.070 


0.303 


0.011 


0.868 


-0.158 


0.020 


0.397 


0.020 


0.868 


m.I.e. 


18 


-0.059 


0.385 


0.016 


0.816 


-0.004 


0.953 


0.718 


0.385 


0.953 


m.m.m 


19 


-0.074 


0.273 


-0.053 


0.435 


-0.134 


0.048 


0.252 


0.048 


0.435 


m.I.e. 


20 


-0.055 


0.415 


-0.022 


0.748 


-0.023 


0.739 


0.634 


0.415 


0.748 


m.I.e. 


21 


-0.058 


0.394 


-0.120 


0.077 


-0.068 


0.314 


0.262 


0.077 


0.394 


m.m.m 


22 


-0.104 


0.126 


0.122 


0.072 


-0.198 


0.004 


0.067 


0.004 


0.126 


m.I.e. 


23 


-0.046 


0.501 


-0.161 


0.018 


-0.135 


0.046 


0.188 


0.018 


0.501 


m.I.e. 


24 


-0.093 


0.173 


-0.118 


0.082 


-0.058 


0.391 


0.215 


0.082 


0.391 


m.I.e. 


25 


0.022 


0.748 


0.071 


0.297 


-0.039 


0.568 


0.538 


0.297 


0.748 


m.I.e. 


26 


-0.031 


0.651 


0.041 


0.544 


-0.210 


0.002 


0.399 


0.002 


0.651 


m.I.e. 


27 


0.006 


0.929 


-0.051 


0.457 


-0.002 


0.976 


0.787 


0.457 


0.976 


m.I.e. 


28 


-0.029 


0.668 


-0.005 


0.938 


0.017 


0.803 


0.803 


0.668 


0.938 


m.m.m 


29 


-0.160 


0.019 


-0.022 


0.743 


-0.020 


0.770 


0.511 


0.019 


0.770 


m.I.e. 


30 


-0 189 


0.005 


0.002 


0.981 


-0.154 


0.024 


0.337 


0.005 


0.981 


m m m 




mean 


mean 


mean 


mean 


mean 


mean 


mean 


mean 


mean 






-0.07884 


0.37010 


-0.05429 


0.44144 


-0.06601 


0.45688 


0.42281 


0.14569 


0.69136 






stdev 


stdev 


stdev 


stdev 


stdev 


stdev 


stdev 


stdev 


stdev 






006784 


0 30264 


0.16079 


0.33413 


0.06379 


0.33785 


0.19490 


0 18499 


0.24739 
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Figure 1. Maximum p-Value for 20 Engines 
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Replications 

Figure 2. Minimum p-Value for 20 Engines 
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Max p-V'alue 



MLE 




Figure 3. Maximum p-Value for 100 Engines 
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Figure 4. Minimum p-Value for 100 Engines 
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V. SUMMARY, CONCLUSIONS AND RECOMMENDATIONS 



A. SUMMARY 

The analytical approach taken to develop this thesis and topics covered has been the 
following: 

1. Develop a program in S-Plus to generate the failure times for an engine with 3 
subsystems, where the distribution of the first failure times of these subsystems is 
Weibull and the engine is not restored to new condition after each engine failure. 

2. Estimate the parameters for the Weibull distribution that fits the generated engine 
failure times with the Maximum Likelihood Estimator and a Modified Method of 
the Moments using commercial spread-sheet software. 

3. Measure the accuracy of these two methods with the Kolmogorov-Smimov Test. 

4. Apply the same analysis to a given set of real engine failure times. 

5. Compare the accuracy of the two methods, and discuss reasons for their 
differences. 



B. CONCLUSIONS 

The two methods in general give good results when they are applied to fit a Weibull 
distribution to the failure times generated in the simulation. The resulting fit for these failure 
times was better with the Maximum Likelihood Estimator more often than with the Modified 
Moment Method. For this reason, and the fact that the MLE method has superior theoretical 
properties in the presence of independence, the author recommends the MLE method to 
estimate the parameters for scenarios similar to those simulated in this Thesis. 

With these conclusions it is important to point out the following considerations; 
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The MMM has the advantage that fewer calculations are necessary to get the results. 
In the case of the Macros developed in this study, only the initial values and two columns of 
calculations are necessary. Four columns of calculations are necessary in the MLE case. 

There was a wide range of values for the parameters of the fitted distributions across all 
the replications. This is displayed in Table 1 1. 



Table 11. Parameters Range 



MLE 

0.961 1<= shape <= 1.5238 
223.28 <= scale <= 330.53 


20 Engines 


MMM 

0.9857 <= shape <= 1.4616 
236.97 <= scale <=329.43 


MLE 


100 Engines 


MMM 


1.1002 <= shape <= 1.3422 




1.1737 <= shape <= 1.4485 


255.64 <= scale <=301.25 




254.07 <= scale <=305.76 



The dependence obtained in the simulation is not as significant as was expected. 



C. RECOMMENDATIONS 

These two methods should be studied when greater dependence exists among the data 
and when the dependency is introduced in a different manner than used in this thesis. 
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APPENDIX A. PROGRAM TO GENERATE FAILURE TIMES 



fiuiction(n) 

{ 

# generate 

# This function generates times randomly in a manner that resembles the way in which 

# engines failures occur in the real world. 

totalreplic <- NULL 
for(count in l:n) { 

# Step 1 

# Three Weibull random numbers are generated each one represents the first failure time 

# that each one of the components could have. These values constitute the vector faults 1. 

# Three Weibull random numbers are generated, each one represents the second failure time 

# that each one of the components could have. These values constitute the vector faults!. 

# Three Weibull random numbers are generated, each one represents the third failure time 

# that each one of the components could have. These values constitute the vector faults!. 

xll <-nveibull(L L2, 900) 
xl2 <- rvveibull(l, 1.3, 800) 
xl3 <- rweibuU(l, 1.4, 600) 
faultsl <- c(xl 1, xl2, xl3) 
x21 <- n\eibuU(l, 1.2, 900) 
x22 <- rv\eibuU(l, 1.3, 800) 
x23 <- rweibullfl, 1.4, 600) 
faults! <- c(x21, x22, x23) 
x31<-rweibuU(l, 1.2, 900) 
x32 <- rv\eibull(l, 1.3, 800) 
x33 <- rv\ eibull(l, 1.4, 600) 
faults! <- c(x3 1, x32, x33) 

failures <- matrL\(c(faultsL faults!, faults!), ncol = 3) 

# Step 2 

# The minimum value of the vector faultsl is selected, this value is the time of the first 

# failure of the engine. The program creates a new vector newfail! that is the difference 

# between the vector faultsl and the minimum value selected, it is clear that this new vector 

# has a component that is 0, this component is replaced in the vector newfail! for the 

# element of the vector faults! that match with the component which had the first failure. 

# Choose the first failure time, that is the minimum 

# time from the veaor faultsl 

fl <- min(faultsl) 

# the next step is to choose the second failure time, to make this first 

# the above value min(faultsl) must be substituted for a new time 

# generated in the vector faults! and precisely in the same component. 

# After this the new minimum is selected 

newfail! <- (faultsl - fl) 
ii <- match(0. newfail!) 
newfail! [ii] <- faults! [ii] 

# Step 3 

# At this point the engine has two components that have accumulated working hours and 

# one that has been just repaired. The vector newfail! has two elements that have the 
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# remainder of the theoretical first failure for these components and one element that is the 
M theoretical second failure for the other component. So, again, the minimum element of 

# tliis vector is selected and it is the second failure time of the engine. 

f2 <- min(newfail2) 
a Step 4 

# The program creates a new vector called ne\vfail3 that is the difference between the vector 

# newfail2 and this minimum value selected, it is clear that this vector has a component that 

# is 0. 

# Choose the third failure time, here there are two possible 

# cases, that the second failure is in the same component that 

# the first one or not, these two condition are programmed in this order 

# 

ne\vfail3 <- (newfail2 - f2) 
iii <- match(0, ne\\faiI3) 
if(ii = iii) { 

ne\\fail3 [iii] <- faults3 [iii] 
f3 <- min(ne\vfail3) 

} 

else if(ii != iii) { 
nevvfail3[iii] <- faults2[iii] 

D <- min(ne\\fail3) 

} 

enginefails <- c(fl, f2, f3) 

#cumenginefails <- cumsum(enginefails) 

#engfails <- matrix(c(failures, enginefails, cumenginefails), b\Tow = 

# T, ncol = 3) 

totalreplic <- rbind(totalreplic, (enginefails)) 

} 

totalreplic 

} 
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Row Names 


VI 


V2 


V3 


1 


88.980852932 


7.374101636 


166.612433434 


2 


57.770945753 


452.164077047 


671 .533220087 


3 


348.363054836 


13.725217444 


33.719418711 


4 


382.257870284 


168.795667771 


119.352260671 


5 


57.145359689 


117.795729216 


228.544288266 


6 


251.488400993 


417.099715694 


206.169150245 


7 


703.336492627 


16.866142294 


168.320579106 


8 


85.174113550 


960.378668893 


561 .250692879 


9 


126.379607253 


166.962891751 


410.881188239 


10 


9.871149363 


30.987527462 


1171.995737937 


11 


184.016586935 


55.183363300 


170.715797827 


12 


19.085994032 


206.177991742 


124.656084035 


13 


171.708369596 


7.411273439 


458.366290751 


14 


957.877021060 


334.941605556 


252.553811069 


15 


88.483061291 


129.532846338 


55.052226186 


16 


97.662694658 


64.764644477 


521.223680534 


17 


213.546020387 


8.955316669 


314.364020817 


18 


420.340591055 


182.036946622 


345.506361499 


19 


46.823610226 


411.857588415 


109.450083602 


20 


25.300667331 


361.149199434 


744.757337329 












Column 1 


Column 2 


Column 3 


Column 1 


1 






Column 2 


-0.099201099 


1 




Column 3 


-0.325682845 


0.177625785 


1 
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Names 


VI 


V2 


V3 


1 


450.634670353 


134.620284266 


763.602013650 


2 


395.264314931 


300.696009190 


201 .433889226 


3 


467.906036969 


301.510351673 


22.166715862 


4 


200.004796884 


261.419709923 


234.01 9837022 


5 


62.187250449 


351.566316602 


293.954533661 


6 


188.420555520 


530.441439163 


446.499347177 


7 


517.045829232 


212.650323103 


4.37229061 2 


8 


104.406621033 


91 .440692508 


286.231002289 


9 


1 85.580658025 


189.658931128 


461.452317203 


10 


278.805459607 


142.687073420 


63.785901768 


11 


242.374076871 


140.003133652 


57.571537504 


12 


191.638384699 


35.982496826 


135.685578893 


13 


115.651429115 


807.515319374 


56.217725189 


14 


318.738490247 


220.511514630 


96.958868093 


15 


389.324028552 


121.699393963 


391 .967785281 


16 


108.549643630 


833.642731237 


587.010089214 


17 


647.283140422 


270.006356536 


600.593877204 


18 


135.277879814 


235.260576261 


294.610909130 


19 


167.565552866 


19.382052477 


80.687152199 


20 


620.417102151 


43.940840941 


348.339229735 


21 


549.652520297 


147.263663645 


402.329146256 


22 


160.858623029 


198.697751449 


131.936073250 


23 


269.748636292 


216.407382425 


272.564254146 


24 


119.882883842 


42.753769287 


568.105938757 


25 


290.367131217 


46.018913125 


470.135356798 


26 


111.526180647 


59.040162467 


442.572648090 


27 


315.869668084 


112.068427109 


799.073849390 


28 


258.473619251 


740.012832046 


730.391096774 


29 


936.074568617 


525.344669063 


53.458836253 


30 


38.275293391 


264.758222487 


109.033459770 


31 


200.399299969 


130.448227050 


564.453719467 


32 


16.387049764 


689.920420726 


124.595174314 


33 


282.346963237 


302.451537677 


17.105846707 


34 


298.005804703 


792.704752916 


277.592370478 


35 


488.070794225 


600.835486551 


77.330402620 


36 


651 .631129729 


186.384979015 


302.942800237 


37 


488.356688854 


282.631614706 


12.094239822 


38 


251.322489919 


639.047414057 


315.553538531 


39 


102.245119335 


121.065856465 


248.600567386 


40 


120.292743494 


281.418670027 


6.512039507 


41 


394.628880256 


329.862437764 


42.172683323 


42 


275.823361976 


400.377074953 


205.279868273 


43 


89.822151642 


568.897212277 


126.945901458 


44 


413.319725633 


280.829982419 


74.621453179 


45 


172.347938766 


138.081515648 


169.690890193 


46 


295.648493372 


49.777202258 


271.146965769 


47 


427.883569557 


120.358845277 


86.140070778 


48 


106.453003164 


336.696295308 


392.945334315 


49 


439.747178352 


145.848786206 


357.670884257 


50 


95.424425759 


158.394437081 


233.214847248 
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Row Names VI V2 V3 



51 


78.318115166 


233.292392917 


150.556841342 


52 


400.322169035 


703.176710554 


131.985760999 


53 


284.082094482 


700.398323845 


54.502213913 


54 


673.813459742 


628.345709323 


27.810511050 


55 


416.157442911 


309.802781935 


61.701311275 


56 


502.727434228 


208.984463635 


7.000471997 


57 


395.249167889 


181.615779362 


246.169600960 


58 


532.055233005 


282.774380172 


51.530633142 


59 


256.526079351 


265.073575146 


366.480188363 


60 


791.013013348 


483.679778572 


383.122947523 


61 


116.674678096 


203.429431816 


501.860086251 


62 


216.019262031 


50.783576381 


346.268316228 


63 


652.055312871 


297.710549841 


19.143684755 


64 


57.143580125 


444.863422850 


53.260846841 


65 


305.430041853 


64.522458661 


344.328881540 


66 


134.875969833 


153.848996778 


74.573546508 


67 


207.795329813 


197.474623424 


36.572283433 


68 


55.079317773 


148.585617029 


153.485311966 


69 


332.356937878 


183.167123218 


2.771444009 


70 


33.292174799 


281.412899135 


521.055842681 


71 


329.237256924 


624.659923834 


410.145804201 


72 


134.999172015 


924.566774196 


959.801158180 


73 


310.283780658 


196.949468983 


61 .420855422 


74 


145.927518733 


98.091275956 


105.358948583 


75 


96.240971550 


448.577731504 


68.349212862 


76 


504.4571 85330 


3.995392340 


54.819812986 


77 


779.185987694 


63.103630524 


510.654006981 


78 


64.295309451 


227.558484641 


681.717084473 


79 


.132.919226467 


80.0901 1 2060 


305.523448635 


80 


225.006666940 


354.405022103 


115.791651829 


81 


322.833702315 


15.741543588 


273.716914308 


82 


552.256221658 


238.880493230 


103.236641646 


83 


234.905588128 


267.774447132 


384.884639091 


84 


241.821803683 


27.660947307 


35.559153707 


85 


360.121368725 


54.400309193 


74.75701 2406 


86 


109.375105806 


169.046509653 


306.549361255 


87 


211.877188712 


54.427577933 


47.413167726 


88 


328.701690509 


233.519929404 


411.204384691 


89 


271 .788277670 


132.777355320 


359.728088427 


90 


54.715071977 


52.466467632 


253.216046834 


91 


157.314280817 


49.880943701 


746.100150202 


92 


93.608619218 


110.699690363 


573.390687287 


93 


218.283430910 


434.282741182 


109.443120744 


94 


1106.554817220 


178.424617471 


107.752314459 


95 


179.703059410 


96.191152623 


63.123472845 


96 


220.753490182 


358.476353401 


438.979104011 


97 


97.987549536 


57.556717482 


152.972644520 


98 


126.992796611 


173.350825362 


1004.747325381 


99 


656.574251587 


250.152057881 


312.989740534 


100 


189.330584986 


142.057479837 


63.629475350 












Column 1 


Column 2 


Column 3 


Column 1 


1 






Column 2 


0.012345873 


1 




Column 3 


-0.151985402 


0 027275574 


1 
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APPENDIX B. MACRO 20 ENGINES 



’ MACRO FOR 20 ENGINES 
' thesis Macro 

’ Macro recorded by FERNANDO GUIMIL 

Sub thesisO 
SheetsAdd 

Sheets("Sheet2").Select 
Range('Al :C2").Select 

AcUveCell FormulaRlCl = "MAXIMUM LIKELIHOOD ESTIMATOR" 

Range("Al:C2").Select 

Rauge("A2"). Activate 

Selection.Font.Bold = True 

With Selection 

HorizontalAlignment = xlCenter 
.Vertical Ahgnment = xlBottom 
.WrapText = False 
.Orientation = xlHorizontal 
End With 
With Selection 

.HorizontalAlignment = xlCenterAcrossSelection 
Vertical Alignment = xlBottom 
.WrapText = True 
.Orientation = xlHorizontal 
End With 

Range("Dl"). Select 

ActiveWindow.SmallScroll ToRight:=10 
With Toolbars(8) 

.Left = 241 
.Top = 228 
End With 

ActiveWindovv.LargeScroll ToRight:=- 1 
Range("Pl:Rl").Select 

AchveCell.FonmilaRlCl = "MODIFIED MOMENT METHOD" 

Range("Pl;Rl").Select 

Range("Ql "). Activate 

Selection.Font.Bold = True 

With Selection 

.HorizontalAlignment = xlCenter 
.VerticalAlignment = xlBottom 
.WrapText = False 
.Orientation = xlHorizontal 
End Witli 
With Selection 

HorizontalAlignment = xlCenterAcrossSelection 
.VerticalAlignment = xlBottom 
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. WrapText = True 
.Orientation = xlHorizontal 
End WiUi 

Range("Nr*).Select 

Acti ve W i ndovv. LargeScroll T oRi ght : =- 1 
With Toolbars(8) 

.Left = 436 
.Top = 226 
End WiUi 

Range("A2‘*). Select 

ActiveCell.FormulaRlCl = "shape parameter" 

With ActiveCell.Charaaers(Start:=l, Length:=15).Font 
.Name = "Arial" 

.FontStyle = "Bold" 

,Size= 10 

.Strikethrough = False 
.Superscript = False 
.Subscript = False 
OutlineFont = False 
Shadow = False 
.Underline = xlNone 
.Colorindex = xl Automatic 
End With 

Range("B3").Select 

ActiveCell.FormulaRlCl = "initial beta" 

Range("B3"). Select 
Selection.Font.Bold = True 
Witli Selection 

.HorizontalAlignment = xlCenter 
. VerticalAlignment = xlBottom 
. WrapTexI = False 
.Orientation = xlHorizontal 
End With 
Range("B5").Select 

ActiveCell.FormulaRlCl = "final beta" 

Range("B5"). Select 
Selection.Font.Bold = True 
With Selection 

.HorizontalAlignment = xlCenter 
.VerticalAlignment = xlBottom 
.WrapTexI = False 
.Orientation = xlHorizontal 
End WiUi 

Range("A7;C7").Select 
ActiveCell.FormulaRlCl = "scale parameter" 
Range("A7:C7"). Select 
Riinge("B7"). Activate 
Selection Font. Bold = True 
With Selection 

.HorizontalAlignment = xlCenter 
.VerticalAlignment = .xlBottom 
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.WrapText = False 
.Orientation = xlHorizonUil 
End With 
Range("A7"). Select 
ActiveCell.FormulaRlCl = "" 

Range("B7"). Select 

ActiveCell.FormulaRlCl = ’’scale parameter" 
Range("B8"). Select 
ActiveCell.FormulaRlCl = ’’landa" 
Range(’’B8").Select 
Seleclion.Font.Bold = Tnie 
With Selection 

.HorizontalAlignment = xlCenter 
.Vertical Alignment = xlBottom 
.WrapText = False 
.Orientation = xlHorizontal 
End With 

ActiveWindow. Small Scroll ToRight:=0 

RangeC’A2:C2’’).Select 

Selection.Copy 

ActiveWindow.SmalJScroll ToRight.-l 1 
Range(’’Q2").Select 
ActiveSheet.Paste 
Application.CutCopyMode = False 
Selection.Clear 

ActiveWindow.LargeScroll ToRight:=-2 

Range(’’A2:C2’’).Select 

Selection.Copy 

ActiveWindow. Small Scroll ToRight:=9 
Range(’’P2").Select 
ActiveSheet.Paste 
Range(’’04"). Select 

ActiveWindow.LargeScroll ToRight:=-l 
Range(’’B3"). Select 
Application.CutCopyMode = False 
Selection.Copy 

ActiveWindow.SmallScroll ToRight:=ll 
Range("Q3"). Select 
ActiveSheet.Paste 

ActiveWindow.LargeScroll ToRight:=-2 
Range('’B5’’).Select 
Application.CutCopyMode = False 
Selection.Copy 

ActiveWindow.SmallScroll ToRight:=12 

Range(’’Q5"). Select 

Application.CutCopyMode = False 

Selection.Copy 

ActiveSheet.Paste 

Range("R7").Select 

ActiveWindow.LargeScroll ToRight:=-2 
Range("B5’’).Select 
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Application. CutCopyMode = False 
Selection.Copy 

ActiveWindow.SmallScroll ToRight;=12 
Range("Q5"). Select 
ActiveSheet.Paste 

ActiveWindovv.LargeScroll ToRight:=-2 
Range("A7:C8").Select 
Application.CutCopyMode = False 
Selection.Copy 

ActiveWindow.SmallScroll ToRight:=l 1 
Range("P7"). Select 
ActiveSheet.Paste 
Range("N7"). Select 

ActiveWindow.LargeScroll ToRight:=-2 

RangeC'Dl 1”). Select 

Application.CutCopyMode = False 

Selection.Copy 

Range("D7").Select 

ActiveSheet.Paste 

Range("B8*’).Select 

ActiveWindow.SmallScroll ToRight:=13 
ActiveWindow.LargeScroll ToRight:=-2 
ActiveWindow.SmallScroll ToRight:=14 
ActiveWindow.LargeScroll ToRight:=-2 
RangeC'El"). Select 
Application.CutCopyMode = False 
ActiveCell.FoimulaRlCl = "Number of observations" 
RangeC'El"). Select 
Selection.Font.Bold = True 
With Selection 

HorizontalAlignment = xlCenter 
. VerticalAlignnient = xlBottom 
. WrapTe.xt = False 
.Orientation = xlHorizontal 
End With 
With Selection 

.HorizontalAlignment = xlCenterAcrossSelection 
.Vertical Alignment = xlBottom 
. WrapTexl = True 
.Orientation = xlHorizontal 
End With 

Rangc("E2"). Select 
ActiveCell.FomuilaRlCl = "n" 

Range("E2"). Select 
Witli Toolbars(8) 

.Left = 536 
.Top = 23 
End With 

Selection.Font.Bold = True 
With Selection 

.HorizontalAlignineiU = xlCenter 
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Vertical Alignment = >JBottoni 
.WrapText = False 
.Orientation = xlHorizontal 
End WiUi 

SheetsC'Sheet 1 "). Select 

Range("B2:B2r').Select 

Selection. Copy 

Sheets("Sheet2"). Select 

Range("D10").Select 

ActiveSheet.Paste 

SheetsC'Sheet 1 "). Select 

Range("C2:C21").Select 

Application.CutCopyMode = False 

Selection.Copy 

Sheets("Sheet2"). Select 

Range("E10").Select 

ActiveSheet.Paste 

Sheets(" Sheet 1 "). Select 

Range("D2:D2 1 "). Select 

Application.CutCopyMode = False 

Selection.Copy 

Sheets("Sheet2"). Select 

Range("F10").Select 

Application.CutCopyMode = False 

Selection.Copy 

SheetsC'Sheet 1 "). Select 

Application.CutCopyMode = False 

Selection.Copy 

Sheets("Sheet2"). Select 

ActiveSheet.Paste 

RangeC'D 10:D29"). Select 

Application.CutCopyMode = False 

Selection.Copy 

RangeC'G 10"). Select 

ActiveSheet.Paste 

Range("E10;E29").Select 

Application. CutCopyMode = False 

Selection.Copy 

Active Window. SmallScroll Down:=16 

Range("G30").Select 

ActiveSheet.Paste 

ActiveWindow.LargeScroll Down;=-l 
Range("F10:F29").Select 
Application.CutCopyMode = False 
Selection.Copy 

ActiveWindow. SmallScroll Down:=35 
Range("G50"). Select 
ActiveSheet.Paste 
Range("D42").Select 
ActiveWindow.LargeScroll Down:=-2 
Range("G10:G69").Select 
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Application. CutCopyMcxle = False 
Selection. Copy 
Range("H10").Select 
ActiveSheet.Paste 
Application.CutCopyMode = False 

Selection.Sort Keyl:=Range("H10"), Orderl:=xlAscending, Header 
xl Guess, OrderCustom:=l, MatchCase:=False, Orientation := _ 
xlTopToBottom 
RangeC'B 10”).Select 
ActiveCell.FormulaRlCl = ’’theta" 

Range("B10").SeIect 
Selection.Font.Bold = True 
With Selection 

HorizontalAlignment = xlCenter 
.Vertical Alignment = xlBottom 
. WrapText = False 
.Orientation = xlHorizontal 
End With 
Range("E3"). Select 
AcdveCelLFormulaRlCl = "60" 

Range("E3").Select 
Selection.Font.Bold = True 
With Selection 

.HorizontalAlignment = xlCenter 
. VerticalAlignment = xlBottom 
.WrapText = False 
.Orientation = xlHorizontal 
End With 

ActiveWindow.ScrollRow= 1 
Range("Gl").Select 

ActiveCell.FormulaRlCl = "Failures times Ti" 

RangeC'Gl"). Select 
With Selection 

.HorizontalAlignment = xlCenterAcrossSelection 
.VerticalAlignment = xlBottom 
.WrapText = True 
.Orientation = xlHorizontal 
End With 

Selection.Font.Bold = True 
Range("G9").Select 
ActiveCell.FormulaRlCl = "Ti" 

Range("G9").Select 
Selection.Font.Bold = True 
With Selection 

.HorizontalAlignment = xlCenter 
.VerticalAlignment = xlBottom 
. WrapTexI = False 
.Orientation = xlHorizontal 
End With 
Selection. Cop>^ 

Range("H9").Select 
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ActiveSheet.Paste 
Range("G9").Select 
Application.CutCopyMode = False 
Selection.Clear 

AcliveW indow. S mall Scrol 1 T oRight :=3 
Range("19"). Select 
ActiveCell.FormulaRlCl = "Ln Ti" 

Rangeri9").Select 
Selection.Font.Bold = True 
Witli Selection 

.HorizontalAlignment = xlCenter 
.Vertical Alignment = xlBottom 
. WrapText = False 
.Orientation = xlHorizontal 
End With 

Range("I10").Select 

ActiveCell.FormulaRlCl = "=LN(RC[-1))" 

Range("1 10"). Select 

Selection.AutoFill Destination:=Range("I10;169"), T>pe 
xlFillDefault 
Range("I10:I69").Select 
Range("H70").Select 
ActiveCell.FormulaRlCl = "SUM" 
Range("H70").Select 
Selection.Font.Bold = True 
With Selection 

.HorizontalAlignment = xlCenter 
.Vertical Alignment = xlBottom 
.WrapText = False 
.Orientation = xlHorizontal 
End With 
Selection. Copy 
Range("170"). Select 
ActiveSheet.Paste 
Range("H71").Select 
Application.CutCopyMode = False 
ActiveCell.FormulaRlCl = "=SUM(R[-61]C:R(-2]C)" 
Range("I71").Select 

ActiveCell.FormulaRlCl = "=SUM(R[-61]C:R[-2]C)" 
Range("172").Select 
ActiveWindow.LargeScroll Do\\ti:=- 3 
Range("J9"). Select 
ActiveCell.FormulaRlCl = "Ti'^B4" 

Range("J9").Select 
Selection.Font.Bold = True 
Witli Selection 

.HorizontalAlignment = xlCenter 
.Vertical Alignment = xlBottom 
.WrapText = False 
Orientation = xlHorizontal 
End WiUi 
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Acti ve Wi ndovv LargeScrol 1 T oRight :=- 1 
Active Wi ndow. Small Scrol 1 ToRight : =2 
RangefJ10").Select 

ActiveCell.FormulaRlCl = "=RC[-2]^R[-6]C[-8]" 

Range(‘'J10").Select 

Selectioii.AutoFill Destinatioa:=Range("J10:J69"), T>pe:=_ 
xlFillDefault 
RangeC'J10:J69’').Select 

Selection.Replace What:="^B*'\ Replacement:="^B4'’, LookAt;=xlPart, 
SearchOrder:=xlByRovvs, MatchCase:=False 
Range("K9'').Select 

Acti ve Window. Smal 1 Scrol 1 ToRight : =2 
ActiveCell.FormulaRlCl = ‘Ti^B6" 

Range("K9").Select 
Selection.Font.Bold = True 
With Selection 

.HorizontalAlignment = .xlCenter 
. VerticalAlignment = xlBottom 
.WrapText = False 
.Orientation = xlHorizontal 
End With 

Range("K 1 0"). Select 

ActiveCell.FormulaRlCl = '‘=RC[-3]^R[-4]C[-9J" 

RangeC'K 10"). Select 

Selection. AutoFill Destination:=Range("K10:K69"), Type:=_ 
XlFillDefault 

RangeC'K 10:K69").Select 

Selection.Replace What:="^B*", Replacement:="^B6", LookAt:=xlPart, 
SearchOrder:=xlByRows, MatchCase:=False 
Range("L9"). Select 

ActiveCell.FormulaRlCl = "LnTI*Ti^B6" 

Range("L9"). Select 
Selection.Font.Bold = True 
With Selection 

.HorizontalAlignment = xlCenterAcrossSelection 
. VerticalAligmnent = xlBottom 
.WrapText = True 
.Orientation = xlHorizontal 
End With 
Range("L9").Select 

ActiveCell.FormulaRlCl = "LnTl * Ti^B6" 

Witli ActiveCell.Cliaracters(Start:=l, Length;=13).Font 
.Nmne = "Arial" 

FontStyle- "Bold" 

Size= 10 

.Striketlirough = False 
.Superscript = False 
.Subscript = False 
.OutlmcFont = False 
.Shadow = False 
.Underline = xlNone 
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Colorlndex = xl Automatic 
End With 

RaiigeC'Ll 0").Select 

AcliveCell.FormulaRlCl = ”=RC[-3]*RC1-1]” 

RangeC'L 1 0 "). Select 

Selection. AutoFill Deslination:=Range("L10:L69"), T>pe 
xlFillDefault 

Range("L10:L69").Select 
Range("M64").Select 
ActiveWindow.LargeScroll Dov\ti:=- 3 
AcliveWindow. SmallScroll Do\vn:= 1 2 
Acli\e Window. ScrollRow = 28 
ActiveWindow.ScrollRow = 1 
AcliveWindow. ScrollRow = 61 
Range("K70").Select 
AcliveCell.FormulaRlCl = "SUM" 

Range("K70").Select 
Selection.Font.Bold = True 
With Selection 

.Horizontal Alignment = xlCenter 
.Vertical Alignment = xlBottom 
. WrapTexl = False 
.Orientation = xlHorizontal 
End With 
Selection.Copy 
Range("L70").Select 
ActiveSheet.Paste 
Range("K71").Select 
Apphcation.CutCopyMode = False 
ActiveCell.FormulaRlCl = "=SUM(R[-61]C:R[-2]C)" 
Range("L71").Select 

ActiveCell.FormulaRlCl = "=SUM(R[-61]C;R[-2]C)" 

Range("L72").Select 

With Toolbars(8) 

.Left = 57 
.Top= 106 
End With 

ActiveWindow.LargeScroll Down:=-4 
Range("M9"). Select 
ActiveCell.FormulaRlCl = "Ti^" 

Range("M9"). Select 
Selection.Font.Bold = True 
With Selection 

.HorizontalAlignment = xlCenter 
.Vertical Alignment = xlBottom 
. WrapTexI = False 
.Orientation = xlHorizontal 
End With 

Range("M 1 0"). Select 

ActiveCell.FormulaRlCl = "=RC[-5]^R[-4]C[4]" 
Range("M10").Select 
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Selection. AutoFill Destination:=Range("M10;M69"), Type:= _ 
xlFillDefault 

RangeC'M10:M69'’).Select 

Selection.Replace What:="^Q*", Replacement:="'X56", LookAt;=xlPart 
SearchOrder:=xlByRows, MatchCase.-False 
ActiveWindow.SmallScroll Down:=8 
ActiveWindow.ScrollRow = 33 
ActiveWindow.ScrollRow = 43 
ActiveWindow.ScrollRow = 48 
ActiveWindow.ScrollRow = 51 
ActiveWindow.ScrollRow = 57 
Range("M70"). Select 
ActiveCell.FormulaRlCl = "SUM” 

Range("M70"). Select 
Selection.Font.Bold = True 
With Selection 

.HorizontalAlignment = xlCenter 
VerticalAlignment = xlBottom 
. WrapText = False 
.Orientation = xlHorizontal 
End With 

RangeCM71").Select 

ActiveCeU.FormulaRlCl = "=SUM(RI-61]C:R[-2]C)" 

Range("M72").Select 

ActiveCell.FormulaRlCl = "mean” 

Range("M73 "). Select 
ActiveCell.FormulaRlCl = "" 

Range("M72").Select 
Selection.Font.Bold = True 
With Selection 

.HorizontalAlignment = xlCenter 
.VerticalAlignment = xlBottom 
.WrapText = False 
.Orientation = xlHorizontal 
End With 

Range("M73 "). Select 
ActiveCell.FormulaRlCl = "=R[-2]C/60" 

Range("M74”).Select 
ActiveWindow.ScrollRow = 1 
ActiveWindow.SmallScroll ToRight:=2 
Range("N9”). Select 

ActiveCell.FormulaRlCl = ”(Ti^6-mean)^2" 

Range("N9”). Select 
Witli Selection 

.HorizontalAlignment = xlCenter AcrossSelection 
.VerticalAlignment = xlBottom 
.WrapText = Tnie 
.Orientation = xlHorizontal 
End WiUi 

Selection.Font.Bold = True 
ActiveWindow.SmallScroll Dov\n:=59 
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ActiveWindovv.LargeScroll Dovvn:=-4 
RangeC'N10’*).Select 

ActiveCeU.FormulaRlCl = "=(RC[-ll-R[63]C[-lir2" 
Range("NlO").Select 

Selection. AutoFill Destination:=Range("N10:N69"), T>pe:= _ 
xlFillDefault 

Range("N10:N69").Select 

Selection.Replace What:="-M*)'^2", Replacemeni:="-M73)^2", LookAt 
:=xlPart SearchOrder:=xlByRows, MatchCase:=False 
Range("09"). Select 

ActiveCeU.FormulaRlCl = "(Ln Ti)^2" 

Range("09").Select 
Selection.Font.Bold = True 
With Selection 

.HorizontalAlignment = xlCenter 
. VerticalAlignment = xlBottom 
.WrapText = False 
.Orientation = xlHorizontal 
End With 

RangeCO10“).Select 
ActiveCelLFormulaRlCl = ••=RC[-6]^2" 

Range("O10").Select 

Selection. AutoFiU Destination;=Range("0 10:069"), T>pe:= _ 
xlFillDefault 

Range("O10:O69").Select 

Range("M70").Select 

Selection. Copy 

Range("N70").Select 

ActiveSheet.Paste 

Range("O70").Select 

ActiveSheet.Paste 

Range("N71").Select 

Application. CutCopyMode = False 

ActiveCeU.FormulaRlCl = "=SUM(R[-61]C:R[-2]C)" 

Range("071").Select 

ActiveCeU.FormulaRlCl = "=SUM(R[-611C:R[-2]C)" 

Range("072").Select 

Active Window. SmallScroll ToRight:=2 

ActiveWindow.LargeScroll Down:=-3 

ActiveWindow.ScrollRow= 1 

ActiveWindow.LargeScroll ToRight:=-l 

With Toolbars(8) 

.Left = 543 
.Top = 37 
End With 
Range("B4").Select 

Active Window. SmallScroU ToRight:=6 
ActiveWindow.SmallScroU Dovvn:=60 
ActiveWindow.LargeScroll Dowti:=-4 
ActiveWindow.LargeScroll ToRight;—! 

Range("B4").Select 
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ActiveCelLFormulaRlCl =__ 

"=((6/(P10"^2))*(R[67]C[13]K(R[67]C[7]^2)/60))*(l/59)r(-0.5)’' 
Range("B4"). Select 
Selection. Copy 

ActiveWindow.SmallScroll ToRight:=10 
Range(”Q4"). Select 
ActiveSheet.Paste 
Application. CutCopyMode = False 
ActiveCelLFormulaRlCl =_ 

•‘=((6/(P10''2))*(R[67]C[-2]-((R[67]C[-8]^2)/60))*(l/59)r(-0.5)'' 
Range("Q5"). Select 

ActiveWindow.LargeScroll ToRight:=-l 
ActiveWindow.SmallScroll ToRight:=-l 
Range("C6 "). S elect 
ActiveCell.FormulaRlCl =_ 
•’=(60/RC[-l])+R[65]C[6]-(60*R[65]C[9]/R[65]C[8])" 
Range("C7‘').Select 

ActiveWindow.LargeScroll ToRight:=-l 
Range(”B9”). Select 

AcUveCell.FormulaRlCl = “=(60/R[62]C[9]r(l/R[-3]C)" 

RangeC'B 1 0"). Select 

Active Window. SmallScroll ToRight:=- 1 

RangeC'Bir‘).Select 

AcUveCell.FormulaRlCl ="=1/R[-2]C" 

RangeC'B 1 2"). Select 
ActiveWindow.SmallScroll ToRight:=l 1 
Range(T6").Select 

ActiveCell.FormulaRlCl = ‘'=(R[65]C[-3])/((R[65]C[-2]r(0.5))“ 
Range("Q9"). Select 

ActiveWindow.LargeScroll ToRight:=-2 
ActiveWindow.SmallScroll ToRight:=9 
RangeC'Q 10*'). Select 
ActiveCell.FormulaRlCl = "tlieta" 

RangeC'Q 10"). Select 
Selection.Font.Bold = True 
With Selection 

HorizontalAJignment = xlCenter 
. Vertical Aligrunent = xlBottom 
.WrapText = False 
.Orientation = xlHorizontal 
Range("Q9").Select 

AcUveCell.FormulaRlCl = "=1/(((1/60)*R[62]C[-4])^(1/R[-3]C))" 
RangeC'Q 1 1"). Select 
AcUveCell.FonnulaRlCl = "=1/R[-2]C" 

RangeC'Q 12"). Select 
End With 



End Sub 



48 



APPENDIX C MACRO 100 ENGINES 



MACRO TO 100 ENGINES 
thesis Macro 

Macro recorded by FERNANDO GUIMIL 



Sub thesisO 
SheetsAdd 

Sheets("Sheet2").Select 

RangeCAl:C2").Select 

ActiveCell.FormulaRlCl = "MAXIMUM LIKELIHOOD ESTIMATOR" 

Range(" A 1 : C2 "). Select 

Range("A2"). Activate 

Selection.Font.Bold = True 

With Selection 

.HorizontalAlignment = xlCenter 
. VerticalAiignment = xlBottom 
. WrapText = False 
.Orientation = xlHorizontal 
End With 
With Selection 

.HorizontalAlignment = xlCenterAcrossSelection 
.VerticalAiignment = xlBottom 
. WrapText = True 
.Orientation = xlHorizontal 
End With 

Range("Dl"). Select 

ActiveWindow.SmallScroll ToRight:=10 
With Toolbars(8) 

.Left = 241 
.Top = 228 
End With 

Active Window. LargeScroll ToRight:=-l 
Range("Pl :Rl").Select 

ActiveCell.FormulaRlCl = "MODIFIED MOMENT METHOD" 

Range("P 1 :R 1 "). Select 

Range("Q 1 "). Activate 

Selection.Font.Bold = True 

With Selection 

HorizontalAlignment = xlCenter 
. VerticalAiignment = xlBottom 
.WrapText = False 
.Orientation = xlHorizontal 
End With 
With Selection 

.HorizontalAlignment = xlCenterAcrossSelection 
.VerticalAiignment = xlBottom 
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.WrapText = True 
.Orientation = xlHorizontaJ 
End With 
RangeC'Nr’). Select 

ActiveWindovv.LargeScroll ToRight:=-l 
With Toolbars(8) 

.Left = 436 
.Top = 226 
End With 
Range(”A2").Select 

ActiveCell.FormulaRlCl = "shape parameter" 

With ActiveCell.Characters(Start:=l, Lenglh:=15).Font 
.Name = "Arial" 

.FontStyle = "Bold" 

.Size = 10 

.Strikethrough = False 
.Superscript = False 
.Subscript = False 
OutlineFont = False 
Shadow = False 
.Underline = xlNone 
.Colorindex = xl Automatic 
End With 
Range("B3").Select 

ActiveCell.FormulaRlCl = "initial beta" 
Range("B3").Select 
Selection.Font.Bold = True 
With Selection 

.HorizontalAlignment = xlCenter 
.Vertical Alignment = xlBottom 
.WrapText = False 
.Orientation = xlHorizontal 
End With 
Range("B5").Select 

ActiveCell.FormulaRlCl = "final beta" 
Range("B5").Select 
Selection.Font.Bold = True 
With Selection 

.HorizontalAlignment = xlCenter 
.Vertical Alignment = xlBottom 
WrapText = False 
.Orientation = xlHorizontal 
End Witli 

Range("A7:C7").Select 
ActiveCell.FomiukiRlCl = "scale parameter" 
Range("A7:C7").Select 
Range("B7").Activate 
Selection.Font.Bold = True 
Witli Selection 

.HorizontalAlignment = xlCenter 
.Vertical Alignment = xlBottom 
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WrapText = False 
.Orientation = xlHorizontal 
End With 

Range('*A7"). Select 
ActiveCell.FomuilaRlCl = "" 
Range("B7").Select 

ActiveCell.FormulaRlCl = "scale parameter" 
Range("B8").Select 
ActiveCell.FormulaRlCl = "landa" 
Range("B8"). Select 
Selection.Font.Bold = True 
With Selection 

.HorizontalAlignment = xlCenter 
. VerticalAlignment = xlBottom 
.WrapText = False 
.Orientation = xlHorizontal 
End With 

ActiveWindow.SmallScroll ToRight:=0 

Range("A2:C2").Select 

Selection. Copy 

ActiveWindow.SmallScroll ToRight:=l 1 
Range("Q2"). Select 
ActiveSheet.Paste 
Application. CutCopyMode = False 
Selection. Clear 

ActiveWindow.LargeScroll ToRight:=-2 
Range(" A2: C2"). Select 
Selection. Copy 

ActiveWindow.SmallScroll ToRight:=9 
Range("P2").Select 
ActiveSheet.Paste 
Range("04").Select 

ActiveWindow.LargeScroll ToRight:=-l 
Range("B3"). Select 
Application. CutCopyMode = False 
Selection.Copy 

ActiveWindow.SmallScroll ToRight:=ll 

Range("Q3").Select 

ActiveSheet.Paste 

ActiveWindow.LargeScroll ToRight:=-2 
Range("B5"). Select 
Application.CutCopyMode = False 
Selection.Copy 

ActiveWindow.SmallScroll ToRight:=12 

Range("Q5").Select 

Application.CutCopyMode = False 

Selection.Copy 

ActiveSheet. Paste 

Range("R7").Select 

ActiveWindow.LargeScroll ToRight:=-2 
Range("B5").Select 
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Application.CutCopyMode = False 
Selection. Copy 

AcliveWindovv. Small Scroll ToRight:=l 2 

Rimge(’'Q5").Select 

AcliveSheel. Paste 

AcliveWindow.LargeScroll ToRight:=-2 
Range(" A7:C8”). Select 
Application.CutCopyMode = False 
Select! on. Copy 

ActiveWindow.SmallScroll ToRight:=l 1 
Range("P7"). Select 
ActiveSheet.Paste 
RangeCN7").Select 

ActiveWindow.LargeScroll ToRight:=-2 
RangeC'Dl l").Select 
Application.CutCopyMode = False 
Selection. Copy 
Range("D7"). Select 
ActiveSheet.Paste 
Range("B8"). Select 

ActiveWindow.SmallScroll ToRight:=13 
ActiveWindow.LargeScroll ToRight:=-2 
ActiveWindow.SmallScroll ToRight:=14 
ActiveWindow.LargeScroll ToRight;=-2 
RangerEl”).Select 
Application.CutCopyMode = False 
ActiveCell.FormulaRlCl = "Number of observations" 
Range("El"). Select 
Selecdon.Font.Bold = True 
With Selection 

.Horizontal Alignment = xl Center 
.Vertical Alignment = xlBottom 
WrapText = False 
.Orientation = xlHorizontal 
End With 
With Selection 

.HorizontalAlignment = xlCenterAcrossSelection 
.Vertical Alignment = xlBottom 
WrapText = True 
Orientation = xlHorizontal 
End With 
Range("E2").Selcct 
ActiveCell.FonnulaRlCl = "n" 

Range("E2"). Select 
Selection.Font.Bold = True 
With Selection 

.HorizontalAlignment = xlCenter 
.Vertical Alignment = xlBottom 
WrapText = False 
.Orientation = xlHorizontal 
End With 



52 



With Toolbars(8) 

Left = 539 
.Top = 30 
End With 

Sheets(" Sheet 1 ”). Select 

Range(”B2:B10r’).Select 

Selection. Copy 

Sheets("Sheet2”).SeIect 

ActiveWindow.SmallScroll Dowti:=1 

Range(”D10”).Select 

ActiveSheet.Paste 

SheetsC’Sheet 1 ”).Select 

RangeCC2:C10r’).Select 

Application. Cut CopyMode = False 

Selection.Copy 

Sheets("Sheet2”).SeIect 

RangeC’E 10”). Select 

ActiveSheet.Paste 

Sheets(” Sheet 1 ”). Select 

Range(”D2 :D 101”). Select 

Application. Cut CopyMode = False 

Selection.Copy 

Sheets(”Sheet2 ”). Select 

Range(”F10").Select 

ActiveSheet.Paste 

ActiveWindow.LargeScroll ToRight:=0 

Range(”D10:D109”).Select 

Application.CutCopyMode = False 

Selection.Copy 

Range(”G10”).Select 

ActiveSheet.Paste 

Range("E10:E109”).Select 

Apphcation. Cut CopyMode = False 

Selection.Copy 

ActiveWindow.SmallScroll Down:=100 
Range(”Gl 10”).Select 
ActiveSheet.Paste 

ActiveWindow.LargeScroll Down:=-6 
Range(”F10:F109”).Select 
AppUcation.CutCopyMode = False 
Selection.Copy 

ActiveWindow.SmallScroll Do\\n:=200 

Range("G210”).Select 

ActiveSheet.Paste 

ActiveWindow.LargeScroll Down:=-l 1 

Range(”G10:G309”).Select 

Application.CutCopyMode = False 

Selection.Copy 

Range(”H10”).Select 

ActiveSheet.Paste 

Application.CutCopyMode = False 
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Selection. Sort Keyl:=Range("H10"), Order l:=xl Ascending, Header 
xlGuess, OrderCustom:=l, MatchCase:=False, Orientation:=_ 
xlTopToBottom 

ActiveWindow. Small Scroll Dowti:=-1 

Range("H2:H6").Select 

ActiveCell.FormulaRlCl = "failures times Ti" 

Range("H2:H6").Select 

Range("H3 "). Activate 

Selection.Font.Bold = True 

With Selection 

.HorizontalAlignment = xlCenter 
.Vertical Alignment = xlBottom 
. WrapText = False 
Orientation = xlHorizontal 
End With 
With Selection 

.HorizontalAlignment = xlCenterAcrossSelection 
. VerticalAlignment = xlBottom 
WrapText = True 
.Orientation = xlHorizontal 
End With 

Range("H9"). Select 
ActiveCell.FormulaRlCl = "Ti" 

Range("H9"). Select 
Selection.Font.Bold = True 
With Selection 

.HorizontalAlignment = xlCenter 
.VerticalAlignment = xlBottom 
.WrapTex! = False 
.Orientation = xlHorizontal 
End With 
Range("19"). Select 
ActiveCell.FormulaRlCl = "Ln Ti" 

Range("I9").Select 
Selection.Font.Bold = True 
With Selection 

HorizontalAlignment = xlCenter 
VerticalAlignment = xlBottom 
.WrapText = False 
.Orientation = xlHorizontal 
End With 

Range("l 10"). Select 

ActiveCell.FormulaRlCl = "=LN(RC|-1])" 

Range("110").Select 

Selection. AutoFill Destination:=Range("1 10:1309"), T}pc:=_ 
xlFillDefault 
Rangc("110:1309").Select 
Witli Toolbars(8) 

.Left = 79 
.Top = 270 
End With 
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Range(“G315”).Select 
ActiveWindovv.SmallScroll ToRight:=2 
Range("I305 *'). Select 
ActiveWindow.LargeScroll Do\vti:=- 1 7 
Range("J9").Select 

ActiveWindow.LargeScroll ToRight:=-l 
ActiveWindovv.SmallScroll ToRight:=2 
Range("J9"). Select 
ActiveCell.FormulaRlCl = ••Ti^B4" 

Range("J9").Select 
Selection.Font.Bold = True 
With Selection 

.HorizontalAlignment = xlCenter 
. VerticalAlignment = xlBottom 
.WrapText = False 
.Orientation = xlHorizontal 
End With 
Selection.Copy 
RangeC’K9“).Select 
ActiveSheet.Paste 

ActiveWindow.LargeScroll ToRight;=-l 
ActiveWindovv.SmallScroll ToRight:=4 
Application.CutCopyMode = False 
AcUveCellFormulaRlCl = 'Ti^B6” 

Range("J10").Select 

AcUveCell.FormulaRlCl = "=(RC[-2]^R[-6]C[-8])" 

Range("J 10*'). Select 

Selection.AutoFill Destination:=Range('*J10:J309"), T>pe:= _ 
xlFUlDefault 
Range(*’J10:J309’*).Select 
Range(*'J300").Select 
ActiveWindow.LargeScroll Down:=-16 
Active Window. SmallScroll Down:=-9 
RangerJ10:J309’*).Select 

Selection. Replace What:="^B*)'*, Replacement:="^B4)", LookAt 
xlPart, SearchOrder:=xlByRovvs, MatchCase:=False 
Range('*K10'*).Select 

AcUveCell.FormulaRlCl = ’'=(RC[-3]^R[-41C[-9])*' 

RangeC'K 1 0"). Select 

SelecUon. AutoFill DesUnaUon:=Range(**K10:K309'*), Type:=_ 
xlFillDefault 

Range(*'K10:K309**).Select 

SelecUon.Replace What:="^B*)", Replacement:="^B6)'*, LookAt 
xlPart, SearchOrder:=xlByRovvs, MatchCase:=False 
Range(*'L8:L9*').Select 
AcUveCell.FormulaRlCl = "Ln Ti * Ti^B6” 
Range('*L8:L9'*).Select 
Range("L9*').AcUvate 
With SelecUon 

.HorizontalAlignment = xlCenterAcrossSelecUon 
.VerticalAlignment = xlBottom 
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WrapText = True 
.Orientation = xlHorizontaJ 
End With 
Range("L9").Select 
Selection.Delete Shift:=xlUp 
Range("L8"). Select 
Selection. Copy 
Range(*'L9"). Select 
ActiveSheet. Paste 
Range("L8"). Select 
Application.CutCopyMode = False 
Selection.Clear 
Range("L9"). Select 
Selection.Font.Bold = True 
With Selection 

.HorizontalAlignment = xlCenter 
. VerticalAlignment = xlBottom 
.WrapText = True 
.Orientation = xlHorizontal 
End With 

RangefLlOO.Select 

ActiveCellFormulaRlCl = ”=RC[-3]*RC1-1]" 
RangeC'LlO"). Select 

Selection. AutoFill Destination:=Range("L10:L309"), Type 
xlFillDefault 

RangeCL10:L309’').Select 
ActiveWindow.LargeScroll Down:=-l 
RangeCM285").Select 
ActiveWindow.LargeScroll Down:=-16 
Active Window. SmallScroll ToRight:=4 
ActiveWindow.LargeScroll ToRight:=-l 
ActiveWindow. SmallScroll ToRight:=8 
ActiveWindow. SmallScroll Down:=33 
ActiveWindow.LargeScroll Down:=-2 
Active Window. ScroURow= 174 
ActiveWindow. ScrollRow = 218 
ActiveWindow. ScrollRow = 252 
ActiveWindow. ScrollRow = 309 
RangeCL310").Select 
ActiveCell.FormulaRlCl = ’’SUM” 

Range("L3 10"). Select 
Selection.Font.Bold = True 
With Selection 

.HorizontalAlignment = xlCenter 
.VerticalAlignment = xlBottom 
.WrapTexl = False 
.Orientation = xlHorizontal 
End WiUi 

Range("L3 1 1"). Select 

ActiveCell.FormulaRlCl = "=SUM(R[-301]C:R[-2|C)" 
Rangc("L3 12"). Select 
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RangerM9"). Select 
ActiveCell.FormulaRlCl =’Ti'X56” 

Range("M9").Select 
Selection.Font.Bold = True 
With Selection 

.Horizontal Alignment = xlCenter 
.Vertical Alignment = xlBottom 
. WrapText = False 
.Orientation = xlHorizontal 
End With 

Range("M10”).Select 

AcliveCell.FormulaRlCl = "=(RC[-5]^R(-41C|4])" 

RangeC'M 1 0").Select 

Selection. AutoFill Destination:=Range("M10:M309”), T>pe:=_ 
xlFillDefault 

RangerM10:M309").Select 

Selection.Replace What:="'X^*)", Replacement:="'X36)", LookAt: 
xlPart, SearchOrder:=xlByRows, MatchCase:=False 
Range("N9”). Select 

ActiveWindovv.SmallScroll Dovvn:=107 
ActiveWindow.ScrollRow= 162 
Active Window. ScrollRow = 214 
ActiveWindow.ScrollRow = 243 
ActiveWindow. ScrollRow = 261 
ActiveWindow.ScrollRow = 282 
ActiveWindow.ScrollRow = 297 
RangefM310").Select 
ActiveCell.FormulaRlCl = "SUM" 

Range('’M3 10’').Select 
Selection.Font.Bold = True 
With Selection 

HorizontalAlignment = xlCenter 
. VerticalAlignment = xlBottom 
. WrapTexI = False 
Orientation = xlHorizontal 
End With 

Range("M3 1 r').Select 

ActiveCell.FormulaRlCl = ’'=SUM(R[-301]C:R[-2]C)” 
Range(”M3 12”).Select 
ActiveCell.FormulaRlCl = "mean" 

Range("M3 12"). Select 
Selection.Font.Bold = Tme 
With Selection 

HorizontalAlignment = xlCenter 
.VerticalAlignment = xlBottom 
.WrapTexI = False 
.Orientation = xlHorizontal 
End With 

Range("M313").Select 
ActiveWindow.LargeScroll Dovvn:=-17 
Range("H3 "). Select 
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ActiveWindow.LargeScroll ToRight:=-l 
Range("E3”).Select 
ActiveCell.FormuiaRlCl = "300" 

Range("E3"). Select 
Selection.Font.Bold = True 
Witli Selection 

.Horizontal Alignment = xlCenter 
.Vertical Alignment = xlBottom 
.WrapText = False 
.Orientation = xlHorizontal 
End With 

ActiveWindow.SmallScroll ToRight:=5 
ActiveWindow.ScrollRow= 144 
ActiveWindow.ScrollRow= 178 
ActiveWindow.ScrollRow = 219 
Active Window. ScrollRow = 264 
ActiveWindow.ScrollRow = 291 
ActiveWindow.ScrollRow = 312 
ActiveWindow.ScrollRow = 274 
ActiveWindow.LargeScroll Dowti:=- 15 
ActiveWindow.SmallScroll Dowti:= 37 
ActiveWindow.ScrollRow = 283 
ActiveWindow.ScrollRow = 312 
ActiveWindow.ScrollRow = 303 
Range("M313").Select 
ActiveCell.FormuiaRlCl = "=R[-2]C/300" 
Range("M314").Select 
ActiveWindow.LargeScroll Down:=-17 
Range("N9"). Select 

ActiveCell.FormuiaRlCl = "(TiX^6-mean)'^2" 
Range("N9"). Select 
Selection.Font.Bold = True 
With Selection 

.HorizontalAlignment = xlCenter 
.Vertical Alignment = xlBottom 
.WrapTexI = False 
.Orientation = .xlHorizontal 
End With 
With Selection 

HorizontalAligmnent = xlCenterAcrossSelection 
.Vertical Alignment = xlBottom 
.WrapTexI = True 
.Orientation = .xlHorizontal 
End Witli 

RangeC'N 1 0"). Select 

ActiveCell.FormuiaRlCl = "=(RC(-21-R[303]C[-1])" 
Range("N 10"). Select 

Selection. AutoFill Destination:=Range("N10:N309"), T>pe: 
xlFillDefault 

Range("N 1 0 : N309"). Select 

Active Wi ndo\\\ S mall Scrol 1 ToRight := 1 
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Range(“N10").Select 

ActiveCell.FormulaRlCl = "=(RC|-l]-R1303]C[-l]r2" 

RangeC'N 1 0"). Select 

Selection. AutoFill Destination:=Range(”N10:N309"), T\pe:=_ 
xlFillDefault 

RangeC'N10:N309").Select 

Selection.Replace Wliat:="-M*)^2", Replacement:- ’-M3 13)^2", LookAt 
:=xlPart SearchOrder:=xlByRows, MatchCase:=False 
Active Window. SmallScroll Down:=-2 
ActiveWindovv. SmallScroll ToRight:=3 
Acti\eWindovv.LargeScroll ToRight:=-l 
Range("B4").Select 

ActiveCell.FormulaRlCl = ”=((6/(PI0^2))*l)" 

Range("B5").Select 

Active Window. SmallScroll Dov\ti:=1 6 
Active Window. ScrollRow= 1 
ActiveWindow.SmallScroll Dovvti:=20 
Active Window. ScrollRow = 612 
Active Window. ScrollRow = 581 
Active Window. ScrollRow = 539 
ActiveWindow. ScrollRow = 452 
ActiveWindow.LargeScroll Down:=-9 
ActiveWindow. ScrollRow = 339 
ActiveWindow. ScrollRow = 309 
Range("H310”).Select 
ActiveCell.FormulaRlCl = "SUM” 

Range("H3 10").Select 
Selection.Font.Bold = True 
With Selection 

.Horizontal Alignment = xlCenter 
Vertical Alignment = xlBottom 
. WrapText = False 
.Orientation = xlHorizontal 
End With 

Range("H3 11"). Select 

ActiveCell.FormulaRlCl = "=SUM(R[-301]C:R[-2]C)" 

Range("H312").Select 

ActiveWindow. ScrollRow = 1 

ActiveWindow. ScrollRow = 465 

ActiveWindow.ScrollRow= 287 

ActiveWindow. ScrollRow = 322 

ActiveWindow.LargeScroll ToRight:=l 

ActiveWindow. ScrollRow = 203 

ActiveWindow. ScrollRow = 269 

ActiveWindow. ScrollRow = 3 

ActiveWindow.SmallScroll ToRight:=-4 

ActiveWindow. ScrollRow = 124 

ActiveWindow. ScrollRow = 223 

ActiveWindow. ScrollRow = 317 

ActiveWindow. ScrollRow = 274 

ActiveWindow. ScrollRow = 276 
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ActiveWindow.SmallScroII Down:=23 

Range("H3 10"). Select 

Selection. Copy 

Range("I3 10"). Select 

ActiveSheet.Paste 

Range("I311").Select 

Application. CutCopyMode = False 

ActiveCeU.FormulaRlCl = "=SUM(R[-301]C:R[-2]C)" 

Range("1312").Select 

ActiveWindow.ScrollRow = 1 

ActiveWindow.SmallScroII Down:=l 

ActiveWindow.SmallScroII ToRight:=2 

Range("09").Select 

ActiveCelLFormulaRlCl = "(Ln Ti)^2" 

Range("09").Select 
Selection.Font.Bold = True 
With Selection 

.Horizontal Alignment = xlCenter 
.Vertical Alignment = xlBottom 
.WrapText = False 
.Orientation = xlHorizontal 
End With 

Range("0 10"). Select 
ActiveCeU.FormulaRlCl = "=RC[-5]'^2" 

Range("0 1 0"). Select 

Selection. AutoFill Destination:=Range("0 10:0309"), T>pe.= 
xlFillDefault 

Range("010:0309").Select 
Range("O310").Selea 
ActiveCeU.FormulaRlCl = "SUM" 

Range("O310").Select 
Selection.Font.Bold = True 
With Selection 

.HorizontalAligmnent = xlCenter 
Vertical Alignment = xlBottom 
. WrapTexI = False 
Orientation = xlHorizontal 
End Witli 

Range("03 1 1 "). Select 

ActiveCeU.FormulaRlCl = "=SUM(R[-301)C:R[-2]C)" 
Range("03 12").Select 
ActiveWindow.SmallScroII ToRight:=-7 
ActiveWindow.LargeScroIl Do\vn:=-15 
ActiveWindow.ScrollRow = 1 
ActiveWindow.SmallScroII ToRight:=l 
Range("B4"). Select 
ActiveCeU.FormulaRlCl =_ 

"=((6/(P10"^2))*(R[3071C[13]-R1307)C[7]/300)/299)^(-0.5) 

Range("B4").Select 

ActiveWindow.SmallScroII ToRight:=6 
ActiveWindow.SmallScroII Do\\n:=-l 
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ActiveCell.FonmilaRlCl =_ 

"=((6/(PI0^2))*(R1307]C|13]-(R[307]C[7]'^2)/3()0)/299r(-0.5) 
Range("B4").Select 
ActiveWindow.ScrollRow = 419 
Active Window. ScrollRow = 300 
ActiveWindow.SmallScroll ToRight;=7 
ActiveWindow.ScrollRow = 1 
Active Window.LargeScroll ToRight:=- 1 
ActiveCell.FormulaRlCl = _ 

"=((6/(PI0^2))*(R[307]C[13]-(R[307]C[7]^2)/300)/299r(-0.5) 
Range("B4"). Select 
ActiveCell.FormulaRlCl = _ 

’•=((6/(P10'^2))*(R[307]C[13]-(R[307]C[7]^2)/300)/299r(-0.5) 
Range("l 10"). Select 
ActiveWindow.SmallScroll ToRight:=2 
Range("O10").Select 
ActiveWindow.LargeScroll ToRight:=-l 
Range("B4").Select 
ActiveCell.FormulaRlCl =_ 

"=((6/(P10^2))*(R[307]C[13]-(R[307]C(7]^2)/300)/299)^(-0.5) 

Range("B4").Select 

ActiveWindow.SmallScroll ToRight:=2 
Range("J10").Select 

ActiveWindow.LargeScroll ToRight:=-l 
ActiveWindow.SmallScroll ToRight:=2 
ActiveWindow.LargeScroll ToRight:=-l 
Range("B4"). Select 

ActiveWindow.SmallScroll ToRight:=2 
Range("JH").Select 

ActiveWindow.LargeScroll ToRight:=-l 
Range("B4"). Select 
ActiveCell.FormulaRlCl =_ 

"=((6/(PI0^2))*(R[307]C[13]-(R[307]C[7]^2)/300)/299)'^(-0.5) 
Range("O10").Select 
ActiveWindow.LargeScroll ToRight:=-l 
ActiveWindow.SmallScroll ToRight:=7 
ActiveCell.FormulaRlCl = "=RC[-6]^2" 

RangeC'0 10"). Select 

Selection. AutoFill Destination:=Range("010:0309"X T>pe:=_ 
xlFillDefault 

Range("O10:O309").Select 
ActiveWindow.LargeScroll Down:=-16 
ActiveWindow.SmallScroll Down:=-12 
ActiveWindow.LargeScroll ToRight:=-l 
Range("B4").Select 

ActiveWindow.SmallScroll ToRight:=3 
ActiveWindow.LargeScroll ToRight:=-l 
ActiveWindow.SmallScroll ToRight:=13 
Range("Q4").Select 

ActiveCell.FormulaRlCl = "=RC[-15]" 

R^ge("Q5").Select 
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AcliveWindow.LargeScroll ToRight:=-2 
With Toolbars(8) 

.Left = 398 
.Top= 139 
End With 

ActiveWindow.LargeScroll ToRight:=l 
ActiveWindow.SmallScroll ToRight:=-9 
ActiveWindow. Small Scroll Do\vn:=3 
Range(”C6").Select 
ActiveCell.FormulaRlCl =_ 

‘'=(300/RC[-l))+R[305]C[6]-(300*R[305]C19)/R[3051C[8]) 
Range("K3 10").Select 
AcdveCell.FormulaRlCl = "SUM” 

Range("K3 10").Select 
Selection.Font.Bold = True 
With Selection 

.HorizontalAlignment = xlCenter 
. Vertical Aligmnent = xlBottom 
.WrapText = False 
.Orientation = xlHorizontal 
End With 

Range("K3ir).Select 

ActiveCell.FormulaRlCl = "=SUM(R(-301)C:R[-2]C)" 

Range("K312").Select 

ActiveWindow.LargeScroll Down:=-17 

ActiveWindow.LargeScroll ToRight:=-l 

Range("C6").Select 

ActiveCell.FormulaRlCl =_ 

"=(300/RC[-l])+R[305]C[6]-(300*R[305]C|9]/R[3051C[8]) 

Range("C6").Select 

ActiveWindow.SmallScroll ToRight:=2 
ActiveWindow. ScrollRow = 432 
ActiveWindow. ScrollRow = 311 
Active Window. ScrollRow = 1 
ActiveWindow. Small Scroll ToRight:=-2 
Range("B10").Select 
ActiveCell.FormulaRlCl = "theta" 

Range("B10").Select 
Selection.Font.Bold = True 
Witli Selection 

.HorizontalAlignment = xlCenter 
-Vertical Alignment = .xlBottom 
WrapTex! = False 
.Orientation = xlHorizontal 
End With 

Range("B9").Select 

ActiveCell.FonnulaRICl = "=(300/R[3021C[9])^(l/R[-3]C)" 

R;mge("B9").Select 

ActiveCell.FormulaRlCl = "" 

Range("C9").Select 
ActiveCell.FonnulaRICl ="" 
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Riinge("B9"). Select 

ActiveCell.FomiulaRlCl = "=(300/R[302JC[9ini/R|-3]C)" 
R;inge("Bl l").Select 
ActiveCell.FomiulaRlCl = "=1/R[-2]C’ 

RaiigeC'B 11"). Select 

Active W i ndovv. Smal 1 Scrol 1 T oRiglit : =9 

Range("Q6"). Select 

With TooIbars(8) 

.Left = 58 
.Top = 25 
End With 
Range("P6"). Select 

ActiveCell.FomiulaRlCl = "((1/299)" 

Range("S3 19"). Select 
Active Window. ScrollRovv = 302 
Range("N3 10"). Select 
ActiveCell.FomiulaRlCl = "SUM" 

Range("N3 10"). Select 
Selection.Font Bold = Tme 
With Selection 

.HorizontalAIignment = xlCenter 
.Vertical Alignment = xlBottom 
. WrapText = False 
.Orientation = .xlHorizontal 
End With 

Range("N311").Select 

ActiveCell.FormulaRlCl = "=SUM(R[-301]C:R[-2]C)" 
Range("N312").SeIect 
ActiveWindow.LargeScroll Down:=-17 
Range("P6"). Select 

ActiveCell.FormulaRlCl = "(M3 1 1)/((N3 1 1)^(0.5))" 
Range("P6").Select 

ActiveWindow.LargeScroll ToRight:=0 
ActiveCell.FormulaRlCl = _ 

"=R[305]C[-3])/(R[305]C[-2])^(0.5)" 

Range("Q9"). Select 

ActiveCell.FormulaRlCl = "=1/(((1/300)*R[302]C[-4])^(1/R[-3]C))" 
Range("Q 1 0"). Select 
ActiveWindow.LargeScroll ToRight:=- 1 
ActiveWindow.SmallScroll ToRight:= 1 3 
ActiveCell.FormulaRlCl = "theta" 

Range("Q 10"). Select 
Selection.Font.Bold = True 
With Selection 

.HorizontalAIignment = xlCenter 
.Vertical Alignment = xlBottom 
.WrapTexI = False 
Orientation = xlHorizontal 
End With 

RangeC'Q 11 "). Select 
ActiveCell.FormulaRlCl = "=1/R[-2]C" 
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RangeC'Q 12'). Select 
ActiveWindovv.LargeScroll ToRight:=-2 
End Sub 
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APPENDIX D. PARAMETER CALCULATIONS FOR THE 20 AND 100 ENGINES 



Row Names V1 V2 V3 



1 


88.980852932 


7.374101636 


166.612433434 


2 


57.770945753 


452.164077047 


671.533220087 


3 


348.363054836 


13.725217444 


33.719418711 


4 


382.257870284 


168.795667771 


119.352260671 


5 


57.145359689 


117.795729216 


228.544288266 


6 


251.488400993 


417.099715694 


206.169150245 


7 


703.336492627 


16.866142294 


168.320579106 


8 


85.174113550 


960.378668893 


561.250692879 


9 


126.379607253 


166.962891751 


410.881188239 


10 


9.871149363 


30.987527462 


1171.995737937 


11 


184.016586935 


55.183363300 


170.715797827 


12 


19.085994032 


206.177991742 


124.656084035 


13 


171.708369596 


7.411273439 


458.366290751 


14 


957.877021060 


334.941605556 


252.553811069 


15 


88.483061291 


129.532846338 


55.052226186 


16 


97.662694658 


64.764644477 


521.223680534 


17 


213.546020387 


8.955316669 


314.364020817 


18 


420.340591055 


182.036946622 


345.506361499 


19 


46.823610226 


411.857588415 


109.450083602 


20 


25.300667331 


361.149199434 


744.757337329 



Column 1 Column 2 Column 3 


Column 1 
Column 2 
Column 3 


1 






-0.099201099 


1 




-0.325682845 


0.177625785 


1 
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A 1 B 1 C 


0 


E 


Ft GJ 


H 


!_J J 






M 1 




0 


P 1 Q 1 R 


1 


MAXIMUM LIKEUHOOD 
ESTIMATOR 




Number 

of 

observati 

ons 


. 


Failures 
times Ti 








K 1 

1 




1 


1 




MODIFIED MOMENT METHOD 


2 


shape parameter 




n 




.i 


















shape parameter 


3 




inrtial beta 
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initial beta 


4 




1 001851 






























1 001851 


5 




final beta 


















“ \ 




j 






1 


final beta 


6 





0 973 


•0 00065 


























7811979 


0.985725: 


7 


scale parameter 

















. 1 


1 








scale parameter 


8 


J 


Ian da 










1 








! 


I 










landa ^ 




9 




0.003973 












Ti 


LnTi 


Ti*B4 


Ti*B€ 


LnTI* ! 
Ti*B6 , 


Ti*Q6 


(Ti* 06- 
mean) *2 


(Ln Ti)*2 




0 00395 




10 




theta 




88 98085 


7 374102 


166 6124 


88 98085 


7 374102 


1.997974 


7 40142 


6 98684 r 


13 95^3 


7 166757 


5141343' 


3 99191 




theta ! 




11 




251 7053 




57 77095 


452 1641 


671.5332 


57 77095 


7 411273 


2.003002 


7 438798 


‘7.021108 


14 063291 


7.202367' 


51397 29! 


4 0120181 




253 1397! 




12 








348 3631 


1372522 


33 71942 


348 3631 


8 955317 


219247 


8 991725 


8 4401^6 


18.503941 


8.6794081 


50729.75 


4 805949) 








13 








3822579 


1687957 


119 3523 


382 2579 


9 871149 


2289616 


9.913067 


9 279395 


21 24625; 


9.553737 i 


50336.66 


5 242343’ 








14 








57 14536 


117 7957 


228 5443 


57.14536 


13 72522 


2.619235 


13.79191 


127881 


33 49504' 


1322152! 


48704 32 


6 860391 








15 








251 48^ 


417.0997 


206 1692 


251 4884 


16.86614 


2825308 


1695456 


' 15 62738 


44 15217 


16 1994'5; 


47398.78 


7 982366 


I 






16 








703 3365 


16 86614 


168 3206 


703 3365 


19 08599 


2.948955 


19 19045 


17 62525 


51 97607 


18 29923; 


46488 89 


8.69633^ 


j 






17 








85 174111 


960 3787 


5612507 


85.1 74 nl 


25 30067 


3230831 


25 4524 


23 18714 


74.91374 


24 16031 


43995 8 


10.438271 


1 






18 








126 3796^ 


166 9629 


410 8812 


1263796 


30 98753 


3.433585 


31 18507 


282439 


9697783 


29.50533! 


41782TV 


11.7895'’ 








19 

20 








9 871149 


30 98753 


1171.996 


9 871149 


33 71942 


3.518074 


3393968 


30 66388 


107.8778 


32 067851 


40741.09' 


12 376841 














184 0166 


55 18336 


170.7158 


184 01661 


46 82361 


3 846388 


47 15812 


42 20483 


162 3361 


44 32199 


35944 4' 


14 7947 
16.06^ 








21 

22 








19065991 


206 1781 


124.6561 


19 08599 


55 0523 


4 008282 


55 46214 


49 40532 


198 0305 


51.99069 


33095.39' 




j 






171 7084 


7.411273 


458,3663 1 


171.7084 


55 18336 


4.010662 


55 5945 


49 51983 


198.6073 


52 11276 


"M050 99t 


16.08541 








23 








957 877 


334 9416 


252 5538 


957 877 


57 14536 


4 045598 


57 57483 


51 23211 


2072645 


53 93868 


32390.42 


16 36686 




1 




24 








8848306 


129 5328 


55 05223" 


88 48306 


57 77095 


4 056486 


58 20629 


51 77774 


210.0357 


! 54 52069 


3218127 


16 45508 




! 




25 








97.66269 


6476464 


521 2237 


97.66269 


64 76464 


4.17076 


65.2665 


57 86709 


241 3497 


61.021291 


rM891.22 


17 395241 








26 

27 








213 546 


8 955317 


314 364 


213546' 


8517411 


4 444698 


85.87764 


75.5421 


335.7618 


79.937^ 


I 23708 04 


1975534' 








1 


1 


420.3406 


182.0369 


345.5064 


420 3406 


88 48306 


4 482811 


89 22022 


78 39614 


351 4351 


82 998261 


22774 98 


20.09561 


1 ' 






28 

29 






- - 


46 82361 
25 30067 


411 8576 
361.1492 


109 4501 
744.7573 


46 82361 
25.30067 


8898085 
97 66269 


4 488421 
4 58152 


89 72309 
98 49432 


7882524 

86.299 


353.8009 

395.3806 


83 45851 
91.47988 


, 22636 27 
20286.93 


2014^ 
20 99032 


1 ! 






lo" 


' 












7.374102 


109 4501 


4 695469 


110 4054 


96 41775 


452 7265 


102.3544 


17307 42 


2204743' 








31 














452 1641 


117.7957 


4768952 


118 84 


103 564 


493 8917 


110.0435 


15343 42' 


22 7429 








32 








. 






13 725^ 


119 3523 


1 4.782079 


120 4133 


104 8953 


5016176 


111.4767 


14990 411 


22 86828 








33 

34 














1687957 


124 6561 


i 4 825559 


125.7744 


109 4281 


528 0518 


116.3583 


13818 88 


23 28602 




















117 7957 


126 3796 


! 4.83929 


1275166 


1109 


536 6772 


117.944 


13448 59 


2341873 


■ 






35 








1 






417 0997 


129 5328 


4 863934 


130.7042 


1135914 


! 552 5011 


120 8442 


i 12784 33 


23 65786 








36 








1 


1 




16 86614 


166 6124 


5.11567 


168 1974 


145.1179 


' 742,3754 


154 8791 


i 6246 208 


26 17008 








37 








[I 






9603787 


166 9629 


5.117772 


168.5518 


145 4149 


! 744 2003 


155 2002 


6195 552 


26 19159 


j- 






38 








! 






166 9629; 168 3206 


5 12587 


169 925 


146 5653 


1 7512749 


156 4442 


6001273 


2627455 








39 








r 






30 98753 


168 7957 


5 128689 


170 4055 


146.9678 


' 753.7523 


156.8794 


5934 025 


26 30345 








40 














' 55 18336 


1707158 


5.14 


172 3475 


148 5943 


n^63.7745 


158 6384 


5666 125 


26 4196 








41 














206 178 


171.7084 


5.145798 


1733515 


149 4348 


' 768.9614 


159 5475 


5530 082 


26 47923 


1 






42 






r 






i 


p7 411273 


182 0369 


^^20421 


183.7987 


1581739 


823.1704 


' 169 0036 


4213102 


27 0838 


- 




1 


43 








i 


1 




334 9416 


184 0166 


t 5215026 


1858013 


159 8474 


833.6083 


i 1708152 


3981 217 


! 27 1965 


1 


1 




44 






r ■ 






! 


129.5328 


206.1692 


; 5.328697 


2082125 


178.5416 


. 951.3939 


1 191.0682 


1835.595 


' 28 39501 








45 

46 




' 










64 76464 
1 8 955317 


206 178 
213.546 


: 5 32874 

; 5,363852 


2082214 

''215.6765 


178 549 
184.7545 


951 4413 
990 9966 


j 191.0763 
i 197.8055 


1834 903 
1303.685 


' 28 39547 
28 77091 








47 


■ 


L — ^ 

i ' 


t 






1 — - 


182 0369 


28 5443 


' 5 43173 


2308534 


197 3685 


1072 053 


211 4932 


502.6057 


! 29 50369 






- — ' 


48 












' 41 1 8576 


251 4884 


5527397 


254 0743 


216 6226 


1197 359 


232 4079 


I 2 262488 


30 55212 








49 












! 


361 1492 


' 252 5538 


' 5.531624 


' 255 1526 


217 5154 


1203 214 


233 3784 


0 284797‘ 


30 59887 








SO 














166 6124 


314364 


; 5.750552 


317.7276 


269 1546 


1547 788 


1 289.5892 


3099 941 


33 06884 








51 














671 5332 


334 9416 


5 813956 


338 5651 


2862824 


1664 433 


t 308 2659 


5528 496 


33 80209 








52 






(- -- - 








33 71942 


345 5064 


5 845011 


349 2642 


295 0648 


1724 657 


J1T7.8483 


7045 298 


34 16415 








53 














119 3523 


348 3631 


5 853245 


352 1573 


297 4383 


1740 979 


320 4 387 


7486 855 


34 26048 
''TT68375 








S4 




r 






1 


L _ 


; 228 5443 


361 1492 


i 5.889291 


365 1 071 


308 0554 


1814 228 


332 029 


9626 936 




1 




55 










i 




206 1692 


3822579 


' 5 946095 


1 386 4877 


325.561 1 


1935 817 


351 1508 13744 92 


" 35 35605 








56 










1 




168 3206 


^410.8812 


' 6.018304 


415 4833 


3492574 


2101 937 


i 377.056 20490.18 


{‘36 21998 






i - — . 


57 














" 561.2507 


i 411.8576 


6 020678 


416 4725 


350 0649 


2107 628 


' 377 93921 20743.81 


36 24856 








58 














' 410 8812 


417.0997 


6 033325 


421 7832 


354 3995 


2138207 


i 382 6805 ' 22132 05 


^40101 




1 




59 








} 






1171996 


420 3406 


6 041065 


425 0666 


357 0786 


2157 135 


[_385 61 1:^3012 67 


^49447 




1 




60 






1- 








1707158 


4521641 


6 114045 


457 3096 


383 3564 


r 2343 858 


1 414 373^ 32566 36 


r37 38155 








61 






j 


1 




1 


124 6M1 


458 3663 


6 127669 


463.5941 


388 4719 
4402143 


1 2380 427 


[ 419 9758 34619 7 


f_37M832 




- 




62 














458 3663 


{521.2237 


6256179 


527 2938 


>754 06 


1 476 6932* 58942 7 


1 3913978 




1 




63 










1 




252 5538 


!_5612507 

|^7r5332 


6336168 
6 509563 


567 8647 


473 0743 


: 2994 639 


512.7587 77755 46 


40 07102 








64 

65 








J 






550523 


t^9 6725 


5632956 


3666 808 611.9437 142908 


* 42 37442 










1 






521 2237 


1 703 3365 


' 6 555835 


^1 9^ 


5892362 


^3862 9351 640 5017 1653151 


42 97898 





66 





A L_B_ 1 _C ] D 1 


r~^ 


1 F 1 G 1 


1 H 1 


1 ' 1 




L... K I 


1 L 


M 


_ N . 1 


1 0 1 P 1 


1 ° 


66 






j 314 364 


1 744 7573: 


i 6 613058 


753.92851 


622 9742 


1 4119.765 


677 6683 


1969196| 


43 73254 




67 






' ^3^5 5064"^ 


^57 877" 


6 864719 


970 1243' 


795 8185 


5463071 


868 4643 


402656.5, 


47 12437 ^ 




6S 






109 4501! 960 3787 


6 867328 


972 6626' 


797 84^ 


! 5479 034 


8707 


405498.9 


'47 16019 




69 






i 7447573 


1171.996 


7 066463' 


1187 424 


968 422 


6843 319 


1 1059 54 


6816622 


49.9349^ 




70 


1 




1 


SUM 


SUM 1 




SUM 


SUM 


SUM 


SUM 


SUM ' 




71 






1 


1 15284 8 


296.1914; 


1 


1 1300826 


1 77584 93 


: 14034.72 


3227644' 


1558.848T 




72 






1 






I 




t * ~ 


mean 








73 






' ' 








' 


! 


^ 233 91211 


' 





67 



Names 


VI 


V2 


V3 


1 


450.634670353 


134 620284266 


763.602013650 


2 


395.264314931 


300.696009190 


201.433889226 


3 


467 906036969 


301.510351673 


22.166715862 


4 


200.004796884 


261.419709923 


234.019837022 


5 


62.187250449 


351.566316602 


293.954533661 


6 


188.420555520 


530.441439163 


446.499347177 


7 


517.045829232 


212.650323103 


4.372290612 


8 


104,406621033 


91 440692508 


286.231002289 


9 


185.580658025 


189.658931128 


461.452317203 


10 


278 805459607 


142 687073420 


63.785901768 


11 


242.374076871 


140.003133652 


57.571537504 


12 


191.638384699 


35.982496826 


135.685578893 


13 


115.651429115 


807.515319374 


56.217725189 


14 


318.738490247 


220.511514630 


96.958868093 


15 


389.324028552 


121.699393963 


391.967785281 


16 


108.549643630 


833.642731237 


587.010089214 


17 


647.283140422 


270.006356536 


600.593877204 


18 


135.277879814 


235.260576261 


294.610909130 


19 


167.565552866 


19.382052477 


80.687152199 


20 


620.417102151 


43.940840941 


348.339229735 


21 


549.652520297 


147.263663645 


402.329146256 


22 


160.858623029 


198.697751449 


131.936073250 


23 


269.748636292 


216.407382425 


272.564254146 


24 


119.882883842 


42.753769287 


568.105938757 


25 


290.367131217 


46.018913125 


470.135356798 


26 


111.526180647 


59.040162467 


442.572648090 


27 


315.869668084 


112.068427109 


799.073849390 


28 


258.473619251 


740.012832046 


730.391096774 


29 


936.074568617 


525.344669063 


53.458836253 


30 


38.275293391 


264.758222487 


109.033459770 


31 


200.399299969 


130.448227050 


564.453719467 


32 


16.387049764 


689 920420726 


124.595174314 


33 


282.346963237 


302.451537677 


17.105846707 


34 


298.005804703 


792.704752916 


277.592370478 


35 


488.070794225 


600.835486551 


77.330402620 


36 


651.631129729 


186 384979015 


302.942800237 


37 


488.356688854 


282.631614706 


12.094239822 


38 


251.322489919 


639.047414057 


315.553538531 


39 


102.245119335 


121.065856465 


248.600567386 


40 


120.292743494 


281.418670027 


6.512039507 


41 


394.628880256 


329.862437764 


42.172683323 


42 


275.823361976 


400.377074953 


205.279868273 


43 


89.822151642 


568.897212277 


126.945901458 


44 


413.319725633 


280.829982419 


74.621453179 


45 


172.347938766 


138.081515648 


169.690890193 


46 


295.648493372 


49.777202258 


271.146965769 


47 


427 883569557 


120 358845277 


86.140070778 


48 


106.453003164 


336.696295308 


392.945334315 


49 


439.747178352 


145.848786206 


357.670884257 


50 


95.424425759 


158,394437081 


233.214847248 



68 



Row Names 


VI 


V2 


V3 


51 


78.318115166 


233.292392917 


150.556841342 


52 


400.322169035 


703.176710554 


131.985760999 


53 


284.082094482 


700.398323845 


54.502213913 


54 


673.813459742 


628.345709323 


27.810511050 


55 


416.157442911 


309.802781935 


61.701311275 


56 


502.727434228 


208.984463635 


7.000471997 


57 


395.249167889 


181.615779362 


246.169600960 


58 


532.055233005 


282.774380172 


51.530633142 


59 


256.526079351 


265.073575146 


366.480188363 


60 


791.013013348 


483.679778572 


383.122947523 


61 


116.674678096 


203.429431816 


501.860086251 


62 


216.019262031 


50.783576381 


346.268316228 


63 


652.055312871 


297.710549841 


19.143684755 


64 


57.143580125 


444.863422850 


53.260846841 


65 


305.430041853 


64.522458661 


344.328881540 


66 


134.875969833 


153.848996778 


74.573546508 


67 


207.795329813 


197.474623424 


36.572283433 


68 


55.079317773 


148.585617029 


153.485311966 


69 


332.356937878 


183.167123218 


2.771444009 


70 


33.292174799 


281.412899135 


521.055842681 


71 


329.237256924 


624.659923834 


410.145804201 


72 


134.999172015 


924.566774196 


959.801158180 


73 


310.283780658 


196.949468983 


61.420855422 


74 


145.927518733 


98.091275956 


105.358948583 


75 


96.240971550 


448.577731504 


68.349212862 


76 


504.457185330 


3.995392340 


54.819812986 


77 


779.185987694 


63.103630524 


510.654006981 


78 


64.295309451 


227.558484641 


681.717084473 


79 


132.919226467 


80.090112060 


305.523448635 


80 


225.006666940 


354.405022103 


115.791651829 


81 


322.833702315 


15.741543588 


273.716914308 


82 


552.256221658 


238.880493230 


103.236641646 


83 


234.905588128 


267.774447132 


384.884639091 


84 


241.821803683 


27.660947307 


35.559153707 


85 


360.121368725 


54.400309193 


74.757012406 


86 


109.375105806 


169.046509653 


306.549361255 


87 


211.877188712 


54.427577933 


47.413167726 


88 


328.701690509 


233.519929404 


411.204384691 


89 


271 .788277670 


132.777355320 


359.728088427 


90 


54.715071977 


52.466467632 


253.216046834 


91 


157.314280817 


49.880943701 


746.100150202 


92 


93.608619218 


110.699690363 


573.390687287 


93 


218.283430910 


434.282741182 


109.443120744 


94 


1106.554817220 


178.424617471 


107.752314459 


95 


179.703059410 


96.191152623 


63.123472845 


96 


220.753490182 


358.476353401 


438.979104011 


97 


97.987549536 


57.556717482 


152.972644520 


98 


126.992796611 


173.350825362 


1004.747325381 


99 


656.574251587 


250.152057881 


312.989740534 


100 


189.330584986 


142.057479837 


63.629475350 












Column 1 


Column 2 


Column 3 


Column 1 


1 






Column 2 


0.012345873 


1 




Column 3 


-0.151985402 


0.027275574 


1 
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B 1 C I 


D 1 


E _| 




G 


H 




J 




»- 1 


M 1 


N 1 


0 1 P 1 Q 1 R 


1 


MAXIMUM UKELIHOOD 
ESTIMATOR 




Number 

of 

observati 

ons 









1 














MODIFIED MOMENT METHOD 


2 


shape parameter 




n 




.J 


failures 
times Ti 








' 








shape parameter 


3 


initial beta 




300 
























initial beta 


4 


1 262411 






























1 26241 




5 


final beta 




















4 


1 








final beta 




6 


1.261028 


-1 9E-05 


. 
























■”17 34948' 


12882321 




7 


scale paranneter 



















_ 




1 




scale parameter 


8 




landa 






























landa 




9 




0 003418 












Ti 


UTi 


Ti*B4 


i 

Ti*B6 


In Ti”^ 
Ti*B6 


ri*Q6 


(Tj-Q6- 1 
mean)*2 


(Ln Ti)*2 




0 003394 




10 




theta 




450.6347^ 


134 6203 


763 602 ’ 


450 6347 ! 2 771444 


1 019368 


3 621402 


3 616305 


3 686347' 


3.717989 


229078T] 


1 039112 




theta 


11 




292 5569' 




3952643 


300 696^ 


201 4339 


3952643 


3 995392 


1 385142 


5746651 


5 735662' 


7944705I 


5 955909 


22840121 


1 9186181 




294 6595 




12 








467 906 


301.5104 


22.16672 


467 906^ 


4 372291 


1 475287. 


6 439285 


6 42617^ 


9 480448 


6 6893171 


2281796' 


2 176472 








13 








200 0048 


261 4197 


234 0198 


200 0048 


651204 


1 8736531 


10.64743 


1061^ 


19 898 


11.175221 


2268263* 


3.510574' 








14 








62.18725 


351 5663 


293.9545 


62 18725 


7000472 


1.9459781 


11 66534 


11.63401 


22 63953 


12.26648^ 


226497^ 


3 786829; 


i . I 






15 








188 4206 


530 4414 


446 4993 


188.4206 


12 09424 


2 492729 


232626 


23.18261 


57.78798 


24 80917 


22273^ 


6 213699 


! 1 






16 








517 0458 


212 6503 


4.372291 


517 0458 


15.74154 


27563031 


32 44627 


32.32293 


89.0918 


34.83971 


2197542! 


7.597208 


f 






17 








104 4066 


91 44069 


286 231 


104 4066 


16.38705 


2 796491 


34 1 3487 


34 00322 


95 08971 


36.69092 


2192057' 


7 820364 








18 








185.5807 


189 6589 


461 4523 


18658071 


17 105851 


2 83942 


36.03582 


35.89471 


101 9202 


38 77717 


2185884 


8 062308 








19 








278 8055 


142 6871 


63.7859 


278 8055 


19 14368' 


2 951973 


41 53768 


41 36859 


122.1189^ 


44 82767 


2168030' 


8.714144 








20 




... “1 




242.3741 


140 0031 


5757154 


242 3741 


19 382051 


2 964348 


42 19167 


42 0192 


124.5595 


45 54801 


2165909 


8.787356 








21 








191 6384 


35 9825 


135.6856 


191 6384 


' 22 16672 


3098592 


49 98356 


49 nool 


154 2169 


54 14712 


2140672 


9.601272 








22 








115.6514 


807 5153 


56 21773 


115.6514 27.66095 


3 320022 


66 104 


65 80144* 


2184622 


72 02099 


“208^89 


11 02254 








'23 




1 




318 7385 


220.5115 


96 95887 


318 7385' 27.81051 


3 325414 


66 55554 


66 25042 


220.3101 


72 52305 


2087238 


11 05838 








24 








389 324 


121 6994 


391.9678 


389 324, 33.29217 


3 505322 


83 52575 


83 12216' 


291 37 


91 43863' 


2032^ 


12 28728 








25 








108.5496 


8336427 


587.0101 


108 5496 35.55915 


3 571198 


9076888 


90.32207 


322.558 


99 53712* 


2009^ 


12.75345 








26 








647 2831 


270.0064 


600.5939! 


647 28311 35.9825 


^583033 


92 13521 


91.68017 


328 493 


101 0663' 


2005578^ 


12.83812 






27 




1 




135 2779 


2352606 


2M,6109^ 


135 2779 36.57228 


3 599291 


94 04576 


9357919 


336 8187 


1032054” 


1999524' 


12.95489^ 






28 








167.5656 


19.38205 


80.68715 


167 5656' 38.27529 


3 644805 


99 60763 


99 10723 


3612265 


109.43751 


1981938 


13.28461 






29 








620 4171 


43 94084 


348 3392 


620 417T] 42.17BF 


' 3.741773 


112.5787 


111.9981 


419.0714 


123.9987 


1941151] 


14 00086 








30 








549 6525 


147.2637 


402 3291 


549.6525 ' 42.7537t1 


r 3 755457 


114.5404 


113 9476 


427.9253 


126204J 


1935011 


14 10346 


i 


1 


! 


31 








160.8586 


1986976 


131 9361 


160 85861 43.94087 


1T782844 


118.5697 


1179516 


446 1925 


130 7361 


19224^ 


^ 30991 


1 


I : 




32 








269.7486 


216 4074 


^72 5643 


269 74861 46.01891 


! 3.829052 


125 6921 


125 0288 


478 7419 


138 7 547 


1900251! 


14 66164 ' 


I 




33 








119 8829 


42 75377 


568 1059 


119 8829 47 41317 


’ 18589 


130 5185 


129 8244 


500 9794 


144 1 9381 


1885285 


14 89111| 






34 








^ 290 3671 


^01891 


470 1354 


290.3671^ 49.7772 


^07557 


138787 


138.0396’ 


539 3977 


153 52131 


1859758 


15.269, 1 


i 




35 




1 




*111 5262 


59 04016! 


442 5726 


111.5262' 49 88094 


3 909639 


139.1522 


138 4025 


541 1038 


1539337 


1858633 


15 28528 1 


I- 




36 

37 




1 


i 


i 315 8697 


112 0684 


799 0738 


315 8697, 50 78358 


[”3 927573 


' 142 3386 


141 5682 


1^56.0193 


157.5314 


1848837 


15 42583 








1 




258 4736 


7400128 


730.391 1 


258 4736, 51 53063 


3 942176 


144 987 


144 1994 


[^8 4593 


160.523' 


1640710 


! 15.54076 






38 








936 0 746 


525 3447 


53 45864 


936 0 746 


52 46647 


3 960174 


1483189 


147.5095] 


r584 16331 


1 1642883’ 


[1830507 


i 15 68298 


I 




39 








38 27529 


264 7582 


109 0335 


' 38 27529" 


53 26085 


1 3 975201 


151 1594 


150 3314 


' 597.5976 


'”T67 4997 


P 1821828 


15.80223 






40 








^ 200 39^ 


4482 


564 4537 


200 3993 


53 45884 


3 978912 


151.8691 


151 0365 


1 600.9608 


168 3022 


M819662 


[J5^174 






41 








16 38705 


'”689 9204' 


n 24 5952 


^638705 


^ 40031 


3 99637 


1552534 


154 3984 


*~617 033 


172 1302” 


1809349 


* 15.97097 






42 








1 282 347 


1 302 4515 


1 17.10585 
277 5924' 


282347 


54 42758 


3.996871 


' 155.3516 


154 496 


^617.5005 
, 618 7606 


M722413 


! 1809050 


] 15.97498 






43 

44 








: 298 0058 


792 7046 


' 298 0058 


54 50221 


3.998241 


155 6206 


154 7632' 


172 5457 


1606232 


! 15 98593 












i 488 0708 


600 8355 


n.3304 


' 488 0708 


54 71507 


4 002139 


156 3882 


155 5258 


1^4358 


1734143 


1805896 


16.01712L 






45 








' 651.6311 


186 385 


' 302 9428^ 


651.6311 54 81981 


^004052 


156 7663 


155 9013 


6242369 


173 842 


1804747 


16.032431 






46 








1 488 3567 


282 6316 


120942"4 


488 356^55 07932 


4 008774 


157 7037 


156 8325 


628 7062 


174 9029 


1801698 


16.070271 


47 








[ 251 3225 


6390474 


315 5535 


251.3225 56.21773 


' 4 029232 


161.8296 


1W9311 


648 4266 


179 5736 


1789380 


M 6 2347^ 






48 








1 102 2451 


121.0659 


246 6006 


^ 1022451^7 14358 


' 4.045567 


1652014 


1642805 


664 6077 


rT83.3925 


1779178 
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APPENDIX E. FUNCTION “tocompare^^ 



fuiiction(v, betamle, ihetamle, betamnim, thetanunni) 

{ # This function calculates the kolmogorov-smimov goodness of tlie fit, 

# and plots of the empical distribution and theoretical distribution 

# for a set of data and tvvo sets of parameters of w eibull distribution. 
win.graphQ 

a <- ks.gof(\\ dist = "weibuir*, shape = betamle, scale = thetamle) 
b <- ks.gof(v, dist = "weibull”, shape = betammm, scale = thetammm) 
par(mffo\v = c(l, 2)) 

cdf compare(v, dist = ’’vveibuir, shape = betamle, scale = thetamle) 
cdf.compare(v, dist = ’'weibull”, shape = betammm, scale = thetammm) 
list(a, b) 

} 
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tocompare(a20. 1 ,0.973,25 1.71 ,0.9857,253 . 14) 

[[ 1 ]]: 

One-sample Kolmogorov-Smimov Test; hypothesized distribution = weibull 
data: v 

ks = 0.0546, p-value = 0.994 

alternative hypothesis. True cdf does not equal the weibull Distn. for at least one sample point. 

[[2]]: 

One-sample Kolmogorov-Smimov Test; hypothesized distribution = weibull 
data; v 

ks = 0.0509, p-value = 0.9977 

Alternative hypothesis: True cdf does not equal the weibull Distn. for at least one sample point. 
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Empirical and Hypothesized weibull CDFs 



Empirical and Hypothesized weibull CDFs 





-solid line h «ni pi'ical d/. vslid line n- Khs am pi'tal d/. 

MLE MMM 

Figure 5. Kolmogorov-Smimov Test; Hypothesized Distribution = Weibull. (a) ks = 
0.0546, p-value = 0.994, (b) ks = 0.0509, p-value = 0.9977 
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tocompare(alOO. 1,1.261,292.56,1.2882,294.66) 

[[ 1 ]]: 

One-sample Kolmogorov-Smimov Test; hypothesized distribution = weibull 
data; v 

ks = 0.03 1 9, p-value = 0.9207 

Alternative hypothesis; Tme cdf does not equal the weibull Distn. for at least one sample point. 

[[ 2 ]]: 

One-sample Kolmogorov-Smimov Test; hypothesized distribution = weibull 
data: v 

ks = 0.0382, p-value = 0.7741 

Alternative hypothesis: True cdf does not equal the weibull Distn. for at least one sample point. 
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Empirical and Hypothesized weibull CDFs Empirical and Hypothesized weibull CDFs 





Figure 6. Kolmogorov-Smirnov Test; Hypothesized Distribution = Weibull. (a) ks = 
0.0319, p-value = 0.9207, (b) ks = 0.0382, p-value = 0.7741 



79 




f M 

■'« '1 






80 



LIST OF REFERENCES 



1 . Caudill, Michel R., “Methods For Performance Goal Setting Of Fielded Jet Engines,” 
Master’s Thesis, Naval Postgraduate School, Monterey, California, June 1995. 

2. Malsbury, John A., “Methods For Determining Goals And Expectations For Fielded 
Jet Engines” Master’s Thesis, Naval Postgraduate School, Monterey, California, 
December 1995. 

3. Mlikan, Mark E., “Methods For Determining Performance Expectations And Optimal 
No Build Times of Fielded Jet Engines” Master’s Thesis, Naval Postgraduate School, 
Monterey, California, June 1996. 

4. Thoman, Darrel R., Lee J. Bain, and Charles E. Antle, “Inferences on the Parameters 
of the Weibull Distribution,” Techtiometrics, Vol. 1 1, No. 3, August 1969, p. 445-460. 

5. Youngreen, “Reliability,” Lecture, Naval Post graduate School, Monterey, California, 
Winter, 1997. 

6. Law, Averill M., and W. David Kelton, Simulation Modeling and Analysis, Second 
Edition, New York: McGraw-Hill, 1991. 



81 




82 



BIBLIOGRAPHY 



Gibbons, Jean Dickinson, Nonparametric Methods for Quantitative Analysis, Second Edition, 
Columbus, Ohio: American Sciences Press, 1985. 

Kececioglu, Dimitri, Reliability and Life Testing Handbook, Vol. I, Englewood Cliffs, New 
Jersey: Prentice-Hall, 1993. 

Leemis, Lawrence M., Reliability: Probabilistic Models and Statistical Methods, Englewood 
Cliffs, New Jersey: Prentice-Hall, 1995. 

Nelson, Wayne, Applied Life Data Analysis, New York: John Wiley & Sons, 1982. 

Statistical Sciences, S-PLUS Guide to Statistical and Mathematical Analysis, Version 3.3, 
Seattle: StatSci, 1995. 

User's Guide - Microsoft Excel, Redmond, Washington: Microsoft Corporation, 1994. 



83 






.4 *|jt| 






^ja 

' VMIV 



A ^•fl'r^'»^'nn^•^»/* |^ »»fv l-'^<4 ff* ' 

-^liA ohH> iJ’ ■ • '■ - > 

. '«.A . j 

I .,<^ <- T-i 



Vi 



■ . 0 



y 

akimii 



iOiy-*ik^ vfmV/ 

* ^ '■ •*'»■•' '' W#vrki 

^1^51 '4 /iLit *i *iUi&ir« j 



‘’srUiH ) ’ Aillfc! u T .>A • 4iwV> 



84 



INITIAL DISTRIBUTION LIST 



No. Copies 

1 . Defense Technical Information Center 2 

8725 John J. Kingman Rd., STE 0944 

Ft. Belvoir,VA 22060-6218 

2. Dudley Knox Library 2 

Naval Postgraduate School 

41 1 Dyer Rd. 

Monterey, CA 93943-5101 

3. Professor W. Max Woods 1 

Code MAAVO 

Naval Postgraduate School 
Monterey, CA 93943-5102 

4. Professor Harold J. Larson 1 

Code OR/LA 

Naval Postgraduate School 
Monterey, CA 93943-5102 

5 ADM DIRECTOR SERVICIOS TECNICOS DE LA ARMADA 1 

ADISTE 
Montalban 2 

28014 Madrid SPAIN 

6. ADM DIRECTOR ENSENANZA DE LA ARMADA 1 

ADEENA 

Montalban 2 

28014 Madrid SPAIN 

7. CAPT DIRECTOR 1 

GIMO 

Arturo Soria 287 

28033 Madrid SPAIN 



85 



8. LCDR Fernando Guimil 2 

GIMO 

Arturo Soria 287 

28033 Madrid SPAIN 



86 



DUDLEY KNOX LiBRARI 
NAVAL POSTGRADUATE SCHOOL 
MONTEREY CA i3943-5101 



DUDLEY KNOX LIBRARY 




3 2768 00338994 1 



